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ABSTRACT 

We present U and B galaxy differential number counts from a field of ^^900 arcmin^, based on 
GOYA Survey imaging of the HST Groth-Westphal strip. Source detection efficiency corrections as 
a function of the object size have been applied. A variation of the half-exposure image method has 
been devised to identify and remove spurious detections. Achieved 50% detection efficiencies are 24.8 
mag in U and 25.5 mag in B in the Vega system. Number count slopes are d\og{N) / dm — 0.50 ±0.02 
for 5=21.0-24.5, and dlog(iV)/dm = 0.48 ± 0.03 for [7=21.0-24.0. Simple number count models are 
presented that simultaneously reproduce the counts over 15 mag in U and B, and over 10 mag in 
Kg, using a A-dominated cosmology and SDSS local luminosity functions. Only by setting a recent 
Zf ~ 1.5 formation redshift for early-type, red galaxies do the models reproduce the change of slope 
observed at Kg — 17.5 in NIR counts. A moderate optical depth [tb — 0.6) for all galaxy types ensures 
that the recent formation for ellipticals does not leave a signature in the U ov B number counts, which 
are featureless at intermediate magnitudes. No ad-hoc disappearing populations are needed to explain 
the counts if number evolution is introduced using an observationally-based z-evolution of the merger 
fraction. 

Subject headings: catalogs — cosmology: observations — galaxies: evolution — galaxies: photometry 
— surveys 
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1. INTRODUCTION 

The first deep CCD measurements and automatic de- 
tection algorithms revealed an excess of faint galaxies 
over the simple extr apolation of the tendency of local 
galaxies ()Tvsonll988j) . Differences between the measured 
surface density of galaxies and the predicted extrapola- 
tion of the local luminosity function (LF) can be related 
to changes in the volume element, to evolution of the 
spectral energy distribution of galaxies, or to the effects 
of merging. Some authors have used LF evolution to 
match number count models to optical data, either in 
density {(j)* ) or in lu minosity {M*) fsee iLillv et aliil991l: 
IMetcalfe et al.ll995l among others), or number evolutio n 
bv collapse l|Glazebrook et a l."1994t iFried et alJl200H) : 
while other works i nsert a pop ulation of blue dwarfs that 
vanishes at z-^0.4 l|Babul fc Reeslll99^ . 

The excess in number counts over non-evolution mod- 
els i s more pronounced a s bluer filters are used (see, 
e.g.. lOdewahn et'aP Il996() . Thus, modeling optical 
and NIR number counts simultaneou sly provides addi- 
tional constraints on galaxy evolution. iBroadhurst et aTl 
(^92) resolved the optical/NIR difference in number 
counts by invoking merging and an enhancement of the 
star formation rate in galaxies at moderate redshifts. 
iGardner et alJ l)1996|) reproduce B, V, I, and K num- 
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ber counts until i? ~ 20 and K 16 mag, using passive- 
evolutionary models with a high B normalization. These 
models, as well as others (e.g., Kauffmann et al."'1994 
Nakata et al . 1999; McCracken et al. 2000; Huang et all 



gta 

2001aj) . had at their disposal shallow or noisy NIR count 



data coming from different sources which often disagree 
with each other. With deeper observations, the discrep- 
ancy between NIR and optical counts became more pro- 
nounced. A degeneracy between the effects of galaxy 
evolution and cosmology gives rise to different inter- 
pretations even of th e same observa t ions. Using data 
from several authors, iPozzetti et alJ l)1996f) found that 
a pure luminosity evolution (PLE) model in an open 
Universe (fi ~ 0) fits number counts, colours and red- 
shi ft distributions reaso nably well in U, hj, rf, I, and 
K. iHuang et al.l l)2001bj) showed that B and K number 
counts from Calar Alto Deep Imaging Survey (CADIS) 
are better reproduced by passive evolution models than 
by no-evolution ones, and that an open Universe is pre- 
ferred to an Einstein-de Sitter (EdS), H. — 1, Universe. 
Totani and collaborators fitted very deep optical and 
NIR data from the Hubble Deep Field (HDF) and the 
Subaru Deep Field (SDF) respectively ijTotani fc Yoshiil 
^000; Totani et al. 2001), including selection effects in 
the model. Although both data sets were well repro- 
duced with a PLE model in a fiat A-dominated Uni- 
verse, optical number counts needed a mild merger rate 
(/i ^ 1), while NIR ones were incompatible with merg- 
ing. Nagashima ct al. (2002) have fitted the same data 
as Totani et al. using a semianalytical model (SAM) 
that includes selection effects. Their results rule out the 
standard CDM, low-density model, and favour a flat A- 
dominated Universe or a low-density, open Universe. 

To some degree, the various interpretations are prob- 
ably affected by field-to-field variations, as well as by 
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TABLE 1 

Depths and Areas Reached by Recent Number Counts Studies in U 

and/or B 



Ref. 


Filter 


Depth 


Area 


Comments 






(mag) 


(arcmin^) 




(1) 


(2) 


(3) 


(4) 


(5) 



Present work 
Capak et al. (20041 



Radovich ct al. (2004) 
Huang et al. (2001bJ 
Kummel fc Wagner" (2001') 
Metcalfe et al. (2001 1 



^^sjjd^^^al^ 

Gardner et al. (20001 



Crawford et al. (20001 
Drive^^aL (1^ 
Hogg et ar7l997T 
Williams et al. (19961 



Metcalfe et al. 




U 24.8 846 2.5m INT/WFC 

B 25.5 888 2.5m INT/WFC 

U 25.4 720 Hawaii HDF-N 

B 26.1 720 Hawaii HDF-N 

U 24.4 2520 VIRMOS Deep Imaging Survey 

B 22.5 1080 CADIS Survey 

Bj 24.5 3857 Northern EcUptic Pole Field 

Buarris 27.5 49 William Herschel Deep Field 

Urgo 26.5 49 William Herschel Deep Field 

F300VK 27.6 5.7 WFPC2/HDF-N 

FAbOW 28.6 5.7 WFPC2/HDF-N 

F300VK 26.9 5.7 WFPC2/HDF-S 

F450VK 28.1 5.7 WFPC2/HDF-S 

u' k g 21.0 1.584 ■ 10^ Sloan Digital Sky Survey 

UV 29.0 ~20 HDFs (STIS/HST) 

FUV 30.0 ~20 HDFs (STIS/HST) 

B 26.0 165.6 2.5m Du Pond Telescope 

B450 29.0 5.7 HST (WFPC2) 

U 25.5 81 5m Hale Telescope 

28.0 5.7 HDF (HST/WFPC2) 

S450 29.0 5.7 HDF (HST/WFPC2) 

B 27.5 19.7 2.5m INT 

28.0 3.5 4.5m WHT 



Note. — Depths have been eonvcrted to magnitudes in the Vega system. 



differing data reduction and analysis techniques. It 
has been a general trait that models that fit opti- 
cal data need to be modified to fit the faint end of 
the NIR counts. One of the key problems revealed 
by recent, high-quality NIR count data is the slope 
change in NIR number counts at K — 17.5 (here- 
after, the "knee"). This feat ure has been re ported 
by several aut hors lIGardner et al. 1993,: Bershadv et al.l 
p98; McCra ckeiiet all 120001: iCristobal- Hornillos et al. 
12003. hereafter CH03D and is now well established. 
While CH03 provide a model that reproduces the Kg = 
17.5 knee, no model has been yet presented that simul- 
taneously reproduces this NIR feature and the counts in 
blue bands, which do not show a knee at intermediate 
magnitudes. 

We are carrying out a deep optical-NIR survey as part 
of a wide project for studying galaxy evolution and for- 
mation, the GOYA Survey^. In this paper we present 
U and B number counts over a ^--^900 arcmin^ area of 
sky, covering one of the GOYA survey fields, the Groth- 
Westphal Strip (GWS). Our U number counts present 
one of the highest product depth x area reached at the 
moment (see Table ^ . The present U and B galaxy 
number counts are complementary to the K5 number 
counts published by our team (CH03), so we have fitted 
a number count model to our optical (U and B) and NIR 
{K) data over the GWS to put new constraints to the 
different ingredients of the galaxy number count models. 

The paper is organized as follows. The GOYA Survey 
and the GWS are described in Comments on the ob- 

4 Known as the 'COSMOS Project' up to 2004 February. GOYA 
Proiect home pag e: 

[http: / /www. iac.es/proyect/GOYAiac/GOYAiac.htmH 



servations are in while reduction is described at 21 
Source extraction and estimation of detection efficiency, 
reliability, and Galactic extinction are presented in 
which summarizes the catalog generation process. In 
^5.41 we explain the procedure to subtract stars counts. 
Final U and B galaxy number counts over the GWS field 
and modeling are presented and discussed in and 
A brief summary is given in SjHl We use a Q,m = 0.3, 
JIa = 0.7, Hq ~ 70 km s^^ Mpc~^ cosmology. 

2. GOYA SURVEY 

The GOY A Survey is described in detail elsewhere (see 
iBalcells eTal. 2002,. CII03, and references therein), so 
we proceed to make a brief introduction of the survey. 
GOYA (Galaxy Origins and Young Assembly) is a wide 
project for studying galaxy formation and evolution with 
EMIR, the NIR multiobject spectrograph t hat will be op- 
erated on the 10 m GTC fsee Balcella 1998tlBalcells et alJ 
[1)00, 2002). 

The GOYA photometric Survey is a multi-color survey 
in six broad band filters {U, B, V, I, J, Kg), cover- 
ing ~0.5 deg^ of sky in several fields, with target depths 
of U=B^V^I^26, and J=i4:=22 (AB mags). Its prin- 
cipal aim is to generate a galaxy database for sample 
selection and characterization for subsequent NIR spec- 
troscopy with EMIR. 

The U and B imaging presented here from INT/WFC 
cover the GWS field (Groth et al. ,1994) . GOYA Survey 
has also reduced and analysed data over this field in NIR 
filters from WHT/INGRID (J and Ks, see CH03), and in 
visi ble filters from HST/W FPC2 {F606W and F8UW, 
see lR,atnatunga et a l."1995l. Originally, GWS field was 
defined as 28 HST/WFPC2 pointings extended along a 
45 arcmin strip, centered at a = 14'^16™38!8 and S ^ 52 
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TABLE 2 
INT/WFC Main Parameters* 



Fig. 1. — Observations of GOYA Survey at the GWS region over 
a DSS image. The 28 HST/WFPC2 pointings defining the GWS_ 
in F6 06W and ^81414^ are shown in grey (see Ratnatunea ct al., 
[T995h . while J and Ks WHT/fNGRID fields are shown in white 
(see CH03 for more information about Kg data). The GWS runs 
diagonally across the wide 45'x45' field marked in black, which 
corresponds to the INT/WFC field (U and B data, presented here). 
Numbers in black indicate positions and orientations of the 4 chips 
of WFC. 



°16'52"(J2000.0) and inclined 40°3'48"to the North. It 
has an area of ~150 arcmin^ of sky. F606W and FSliW 
data in the field are provided by the DEEP database^ 
(sec PhiUiDS ct al. 1997; Simard et al. 2002; Weiner et all 
120051) . as well as morphology and photometry. Exposure 
times were 4,400 s in F8UW and 2,800 s in i^606W^ for 
27 pointings, and 25.2 ks in both WFPC2 filters for a sin- 
gle pointing. In Figure ^ covered areas in the available 
six filters of the GOYA survey are plotted over a DSS® 
image of the GWS sky region. Compared to other exist- 
ing optical-NIR surveys, GOYA offers a notable increase 
in the depth X area product in several filters, compiling 
complementary photometry in six optical-NIR bands, 
and morphological and surface brightness information 
from high-resolution IIST/WFPC2 images. 

3. OBSERVATIONS 

U and B observations were obtained during one run 
in May 2002, using the Wide Field Camera (WFC) 
mounted on the prime focus of the 2.5m Isaac Newton 
Telescope (INT) at Roque de Los Muchachos Observa- 
tory, in La Palma. The camera consists of a mosaic of 
four CCDs (see Figure[lJ, each of them with 2,048x4,096 
pixels, giving an irregular field of view of approximately 
34'x34', and a pixel scale of 0.3334'7pixel. The main 
camera and filter characteristics are listed in Table [3 
The average interchip spacing is ~ 1". In order to cover 
these gaps and to facilitate bad pixel correction in the 
final stacked images, we dithered the exposures by ~10" 
in each band, but only in the N-S direction. Therefore, 

^ DEEP Project Home Page: 
[http: / /deep. ucolick.org/ 

^ Dig itized Sky Survey ( DSS): 

[http: //a rchive. stsci.edu/ds s/index.htrnl| 



Parameter 


Value 






„ ,, . 

C_iOiiGCtrlIlg 3,rG3, 


5.07 m2 






Focal ratio at WFC 


3.29 






Field of view 


Irregular, ~34'x34' 






Detector type 


4 EEV CCDs 






Detector format 


2048x4096 pixels 






Pixel size 


13.5xf3.5 fj.ra 






Pixel scale 


0.3334 arcsec pixel" ^ 






Field co-planity 


±20 fx 






Readout time 


56 s (full camera) slow mode 




Cosmic ray counts 


~ 2 000 per hour per chip 




(at sea level) 








Chip planity 


6-10 fim 






Bad pixels 


< 2% 






Parameter 


CCD#1 CCD#2 


CCD#3 


CCD#4 


■ / — A T\TT — 1 \ 

Gam (e ADU '-) 


2.8 3.0 


2.5 


2.9 


Bias (ADU) 


f527 1590 


1623 


1644 


Readout noise (e~ ) 


6.4 6.9 


5.5 


5.8 


Quantum efficiency 


67% 72% 


62% 


61% 


{U, @ -f20°C) 








Quantum efficiency 


80% 87% 


80% 


78% 


(B, @ -120°C) 








Dark current 


3.8 3.3 


2.9 


2.0 


(ADU/hr, @ -I20°C) 








Filters 


Peak (A) 


Width (A) 


Urgo 


3518 


638 




Bkpno 


4407 


1022 





INT/WFC information: 

|http:/ /www. ing.iac.es/Astronomy /instruments/ wfc/index.htm]| 



gaps between CCDs ^1, #3 and were removed, but 
not the spacing between CCD^2 and the others (see Fig- 
ure^. The dithering was set to a low value to maximize 
the area of maximum exposure time and to asses this 
last to be uniform. Each image was exposed for 1,800 
s, except one B exposure of 1,300 s, for a total integra- 
tion time of 14,400 s in [/ and 10,300 in B. Effective 
exposure times must be lower due to the loss of light be- 
cause of high cirrus at the beginning of the night. Dome 
and twilight flat-field images were obtained for mapping 
different pixel responses, and zero exposure frames were 
taken to estimate the bias structure in each CCD. 

The North-East corner of the held of view suffers from 
serious vignetting. We therefore offset our pointing, and, 
as a result, miss a fraction of the WFPC2 frame corre- 
sponding to the South- West end of the GWS. 

Atmospheric turbulence produced an average PSF of 
FWHM-1.3" in U and -1.2" in B in the final stacked 
images, stable over the whole run, with small ffuctua- 
tions across the field. The maximum ratio between the 
PSF FWHMs of the inner and outer parts of the mosaic 
was ^ 1.2 in both bands. Standard photometric star 
fields (Landoh 199^ were observed during the night at 
different airmasses in order to correct final images for at- 
mospheric extinction and to determine the photometric 
zero point for each filter. Attained limiting magnitudes 
at 50% detection efficiency are 24.8 mag in U and 25.5 
mag in B in the Vega system (see f|5lfor a description of 
the efficiency and reliability analysis). 
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4. REDUCTION 
4.1. Pre-reduction 

Basic reduction was carried out using a package spe- 
cially designed by us for reducing INT/WFC images, and 
based on t he MSCRED ^ package in IRAF (Valdcs 2002, 
l2nnTl rmflfll fVa,ldes Todvl flflOSV CCD mosaic 

exposures were bias- and dark-corrected using overscan 
columns, because the bias structures were lower than 
0.1%. The ING group® have reported departures from 
hnearity of ~2% starting at ~50,000 ADUs. CASU INT 
linearity coefficients^ were used to correct for this, which 
are known to be very stable over periods of years, and are 
precise to 0.2% in the 0-50 K count range. Images were 
corrected for pixel-to pixel-response using flat-fields from 
the combination of twilight and dome exposures. Vignet- 
ing was completely corrected in CCD^2, but residuals 
at the North-East corner of CCD#3 remained after flat- 
fielding (-5% the sky level in B band and ~20% in U). 

After flat-fielding, a super-flat image could not be con- 
structed from our data because of scattered light from 
internal optics and saturated stars, which introduced dif- 
fuse variable patterns of more than 50% of the sky level 
over scales of several arcmins. The 10" of our dither 
pattern was insufficient for filtering out such diffuse pat- 
terns when constructing the superflat. Diffuse light pat- 
terns in w ide field cameras may affect the photometry. 
iLauer &^a ldcs (1997) found that diffuse light affects on 
small scales when combining images, and Capaccioli et al. 
(2001) 1° found that it inserted a ~3% error in the wide 
bands of the Capodimonte Deep Field, [jvlanfroi d et al.l 
(|2002t) showed that the zero point changes across the 
WFI mosaic were a consequence of scattered light, in- 
serting an additional error of .1 mag in the ?7 zero p oint 
of the VIRMOS photometry l|R,adovich et al J 120041) . In 
addition to the above problems, we argue that variable 
sky levels insert errors in the structural parameters es- 
timated from isophotal analysis. Thus, removing diffuse 
light in mosaics i s need ed in order to get a reliable pho- 
tometry. IValdesI ()2000j) argued that a theoretical model 
of the response of the detector-telescope system was nec- 
essary to remove diffuse light in the NOAO mosaic. We 
built a model surface of the sky by fitting 1-D splines 
to image in both directions consecutively, using rejection 
algorithms to suppress objects from the fitting. A first 
aproximation to the diffuse light pattern was obtained by 
fitting 1-D splines to all the rows of each image, indepen- 
dently from row to row. The posible row-to-row residuals 
of this first model were smoothed out when we performed 
a second fit to all the columns of this first model. The 
result was the desired model surface of the sky, with low 
RMS in smaU areas of 5 x 5 pixels (RMS< 0.01). After 
substracting each sky model to its corresponding image, 
diffuse light residuals were reduced to less than 1% of 

MSCRED Home Page: 
|http: //iraf .noao.edu/ sc ripts /irafref ?mscred| 

" INT Home Page: 
IStp: / /www. ing.iac.es 

^ CASU INT Wide Field Survey Home Page: 
[http: / /www. ast.cam.ac.uk/ ^wfcsur/index.php 

^" Capaccioli, M., et al.2001, The Capodimonte Deep Field: 
Data reduction and first results on galaxy clusters identifica- 
tion, Osservatorio astronomico di Capodimonte (Napoli: Istituto 
Nazionale di Astrofisica) , _ 
|http:// www.na.astro.it/oacdf/OACDFPAP /OACDFPAP.htmI| 



the sky in the whole field of all the exposures (see Figure 
EJ. We verified that aperture photometry of stars was 
unaffected to better than 0.007 mag typically, obtaining 
higher differences in the regions with the most complex 
scattered light structures (< 0.015 mag). 

From regions of images which were not affected by vi- 
gneting or diffuse patterns, it could be deduced that the 
illumination pattern was less than 0.5% the sky level. 
Diffuse pattern residuals are of the same order, so we did 
not correct for it. 

4.2. Photometry 

Several Landolt (1992) fields were taken for photomet- 
ric calibration. Colors of the Landolt stars covered a wide 
range (—1.2 < U — B < 2.0). Fitted calibration equa- 
tions included zero point, atmospheric extinction, and 
color terms, as follows: 

mv = U + ua + ui- X + U2 - (U - B) (1) 
mB=B + bo + bi-X + b2-{U-B), (2) 

where mu and rriB are the instrumental magnitudes in 
the U- and i?-bands respectively; U and B are the John- 
son magnitudes; uq and bg represent the zero points; ui 
and bi the extinction coefficients; U2 and &2 are the color- 
term coefficients; and X represents the airmass. 

Standard stars were positioned in the center of the 
WFC field, on CCD#4, which is free of vigneting. This 
calibration applies directly to the rest of CCDs because 
all the CCDs have been converted to "mean count" units, 
multiplying the flat-fields by the constant factor N/{G) 
; where N is the mean level of the flat-field and (G) is 
the mean gain of the 4 CCDs (see their values in Ta- 
ble Our photometric solution was derived from star 
flelds exposed during the second half of the night, given 
that high cirrus was present during the first half. It was 
applied to the entire dataset by previously scaling sci- 
ence frames from the first half of the night to a reference 
exposure from the second half. Photometric calibration 
results are shown at Table El Final RMS residuals are 
0.09 mag for the f7-band and 0.06 mag for the i?-band. 
Color terms and zero points were similar to those from 
the web page of the INT Wide Field Survey. Estimates 
for the U and B extinction terms derive d from the the- 
oretical extinction curve for La Palma (f Kind IT985I) . as 
well as the extinction measurements of the Carlsberg^^ 
and Mercator^^ telescopes for the night of the run, were 
similar to those obtained from our fits. 

We applied zero-point and extinction correction to all 
of our sources, but color terms were only applied to those 
sources which had counterparts in both fllters. Photo- 
metric errors were obtained by quadrature-sum of photo- 
metric calibration error, Poisson noise and Galactic ex- 
tinction errors. Error values for sources brighter than 
the 50% limiting magnitude are typically '--^O.IO mag in 
U and ~0.05 mag in B. 

4.3. Astrometry 

Astrometry was performed using IRAF astrometric 
tasks. The Guide Star Catalog II" (GSC-II) provided 

Carlsberg Telescope Web: 
http:/ /www. ast.cam.ac.uk/ ^dwe /SFR/camc-extinction.html] 
Mercator telescope Web: ^^^^^^^^^^^^^^^^^^ 

http: / / www. mercator .iac.es/ extin ction/extin_previous.html| 

G uide Star Catalog II: 
http: / / www-gsss.stsci.edu/gsc / gsc2 / GSC2home.htm| 
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Fig. 2. — Left: Typical diffuse structure in one of the B exposures, due to saturated stars' light and reflections in the internal optics 
of the telescope. Beginning from the left, CCDs ID numbers are: CCD#:3, CCD#:4, CCD#:1, and CCD#2. Vigneting at left corner of 
CCD#:3 could not be completely removed by flat-fielding. Right: The same B pointing once the diffuse pattern was subtracted. Grey 
scaling in both figures is not the same, in order to emphasize small residual structures near very bright stars' haloes that resulted after 
diffuse pattern removal, of ~1% sky level. 



TABLE 3 

Photometric Calibration Coefficients for U and B 



Filter Zero Point'' Extinction*^ Color Term'^ RMS'^ X"= 

U 23.64±0.16 0.52±0.11 -0.13±0.04 0.093 1.24 
B 25.34±0.12 0.25±0.09 -0.09±0.03 0.062 1.62 

^Magnitudes per ADU s~^, in the Vega System 

Extinction coefhcient (mag airmass"^) 

Color coefficient (adimensional) 
'^The rms of photometric calibration fit 

'^Airmass of the combined images taken as reference in each 
fiher 



the coordinates of stars in the field. As with most wide- 
field devices, the WFC is known to have an important 
"pincushion" distortion that scales as r^, where r is the 
distance to the optical axis llTavloill200n() . The World 
Coordinate System (WCS) inserted by the telescope in 
the image headers included linear terms only, with a code 
which was unrecognizable by IRAF tasks. Therefore, af- 
ter creating a new primary linear WCS in each CCD, we 
proceed ed to describe WFC distortions using msctpeak 
i|Valdesl 1200(1 ITqQTI) . with the TNX projection, which 
includes high-order polynomial terms to the tangential 
projection fit. The msctpeak task has the disadvan- 
tage of fitting distortion te rms in the transformed co- 
ordinate plane, y and f (see iCalabretta fc Greiseijl2003 [ 
iGreisen &: Calabrettall2003 for a detailed description of 
WCS coding), where the dependence is diluted in a 
complex combination of cross-terms. As the fitting algo- 
rithm is not able to give the right weight to each term, 
astrometry at the CCD edges is not improved by the in- 
clusion of higher-order terms. This fact became evident 
when we combined the 4 separate WFC frames into sin- 
gle images: residuals of up to '--^l" changed sign abruptly 
at the chip-to-chip transitions. Our conclusion was that 
a single fit to the entire field was needed. Combined 
single images were created using the relative positions 
and rotations between the four WFC CCDs reported in 
Uavlor (2000). Discontinuities in the astrometric solu- 
tions indicated that error s of ^1" are pr esent in the chip 
separations presented by iTavloil l)200Cl() . We estimated 
corrections by procedure, offseting the CCD relative po- 



sitions until we found a combination which minimized the 
astrometric RMS in the whole field. We list corrected val- 
ues in Appendix^ The achieved RMS astrometric error 
is down to 0.3" over the entire ^36'x36' field. We have 
also detected a global rotation of our WCS with respect 
to that which is in the DEEP database for the F606W 
and F814W HST/WFPC2 images of the GWS, which is 
not related with the different celestial reference frames. 
Our coordinate system is rotated ~ 0.075° in the tan- 
gential plane to the sky respect to the DEEP coordinate 
system, to the North-East direction. This rotation is cen- 
tered at a = 14.28333'^, S = 52.35°. The problem could 
be due to the small number of stars that were used for 
calibrating the HST images (just 4 for the entire GWS, 
Groth priv. com.). Our catalogues include source coor- 
dinates in both the DEEP and GSC-II systems, in order 
to facilitate cross-correlation with the DEEP database 
information. 

4.4. Coaddition 

Once the astrometric solution and the flux scaling fac- 
tors were computed, MSCRED tasks were used for stacking 
images. After subtracting the background, images were 
re-pixeled in order to create single images of the expo- 
sures on an uniform grid on the sky, free from distortions. 
Resampling was done using a "sine" function multiplied 
by an interpolation kernel, which maintains the statisti- 
cal characteristics of the sky noise (Valdcs 1998). After 
removing atmospheric refraction effects and applying the 
flux scaling factor (see tl4.2() , stacking of images in U and 
B was done by computing the exptime-weigthed aver- 
age of the images, a process which preserves the Pois- 
sonian nature of the sky noise (Valdcs 1998). MSCRED 
algorithms were used in order to reject cosmic rays and 
satellite tracks, masks of saturated defective pixels, and 
bleed trails. In Figure|31we show final J7-band combined 
image, and its exposure time map. The B stacked image 
is similar, but ^0.7 mag deeper. Both images cover an 
irregular-shaped field over the GWS of ~40'x40', with 
an inter-ga,p non-covered zone between CCD#2 and the 
other three. The main U and B image characteristics are 
fisted in Table H 

5. SOURCE EXTRACTION 

A first estimate of the detection limiting magnitude at 
3cr level for point-like sources can be done considering 
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TABLE 4 

Main Characteristics of the Final Stacked U and B Images 



Filter 


Dithering 


Tot. exp. 


FWHM 


Limiting magnitude 


re 


Det. Thresh. 


Area for counting 








(arcsec) 


(mag) 


(arcsec) 


(sky a) 


(arcmin^) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


U 


8 X 1800 


14 400 


1.3 


24.8 


1.5, 3 


0.6 


846 


B 


1300 + 4 X 1800 


10 300 


1.2 


25.5 


1.5, 3 


0.6 


888 



Note. — Col. (6) shows the effective radius, rc, for group division to define "point", "intermediate-size", and "large" objects 
(a detailed description is at i|5. II and i|5.2l . Col. (7) is the detection threshold used in SExtractor in sky a units. The final 
maximum-exposed area over which number counts have been computed in each band is listed in Col. (8) (see JbI . 




Fig. 3. — Upper panel: Final U stacked image over the GWS 
field. B final image is similar. Lower panel: U exposure-time map 
produced after stacking. 



that: 



mag(3CT) = flo — 2.5 • log{3a^/A), 



(3) 



where ao is the zero point in each band, and A is the 
area of the selected aperture. We obtained t/um ^24.4 



and Slim ~25.9 at 3tT for an aperture of 1.0". But a more 
accurate determination of the photometric limits of the 
survey must be done. 

Catalogues we r e ob tained using SExtractor 
(jBertin fc Arnoutl I1996() . which basically considers 
as a detection every group of connected pixels above a 
fixed detection threshold, DETECT_THRESH, after filtering 
the image with a detection kernel. CH03 indicate that 
the detection kernel and the minimum area are well 
constrained as a function of seeing. The minimum area 
was fixed to the area of a symmetrical source with a 
diameter equal to one fourth the FWHM of the stars 
in the image. On the other hand, fixing the detection 
threshold is more difficult because number of detected 
sources will increase as we lower the threshold, but so 
will the spurious detections (noise peaks and source 
substructures identified as sources). A compromise 
betwen maximization of detections and minimization of 
the spurious fraction must be found. 

Incompleteness effects are quantified through Monte 
Carlo methods: the behaviour of the detection algo- 
rithm is characterized studying how it detects an a priori 
known source distribution, inserted in known positions in 
the science frames, or in synthetic images that simulate 
the noise characteristics of the science ones. Previous 
studies of this nature find that efficiency and reliabil- 
ity estimates based on synthetic images tend to over- 
estimate the efficiency and underestimate the spurious 
fraction, probably due to overly regular shapes of artifi- 
cial sources, and to the fact that the real sky noise is not 
extrictly Poissonian (Bershady et al. 1998; CH03). Nev- 
ertheless, the study with synthetic images is useful to 
narrow down the range of DETECT_THRESH values for the 
more detailed study with the science images. Thus, we 
have measured efficiencies and spurious fractions first on 
synthetic images (see H5.1|l . and after that, we have car- 
ried out a more detailed study of incompleteness effects 
using the science frames, which is going to be taken as 
the definitive one for completeness correction ( ti5.2|) . For 
spurious characterization when using the science images, 
we have slightly modified the method of CH03 ( ^5.2.2^ . 

SExtractor is more efficient at detecting high surface 
brightness sources, i.e., compact sources have a higher 
probability of being found than extended sources at the 
same magnitude. We determine this dependence in the 
definitive study (with science frames) by carrying out 
the efficiency analysis over three different groups of ob- 
ject sizes, as CH03 describe. Bins are defined using the 
half-light radius of stars in the images, re, stars: objects 
with To < 1.5 To, stars are "point-like" sources, those with 
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1-5 Te, stars < f'c < 3 To, stars are " intermediate-sized" ob- 
jects, and those with re > 3 rc, stars are "large" objects. 
Therefore, our efficiency study using science frames mea- 
sures detection efficiency not only as a function of the 
source magnitude, but also of its size f H5.2.1|l . 

Differential number counts are usually corrected for 
efficiency using the efficiency matrix, P{m\-n \mrmt), or 
the e fficiency function, E{mout)- Following lYan et al.l 
(jl998j) . the element P(TOin; Wout) = ^'in.out of the effi- 
ciency matrix is defined as the probability of objects with 
an original magnitude mjn to be detected with magnitude 
n^out: 

_ iVdct(™in,"lout) 



P(min; mout) 



A'"orig(TOin) 



(4) 



where -/V(jct(™in7 "^out) is the number of objects with an 
original magnitude mjn which are recovered with magni- 
tude mout, and A^orig(™in) is the total number of objects 
that originally had magnitude mm- Notice that this def- 
inition is independent of the original source distribution. 
This is due to the fact that the probability for an ob- 
ject with an original magnitude to be recovered by 
SExtractor in mout only depends on its size, its origi- 
nal input magnitude min, and the noise characteristics of 
the image, but not on the total number of objects with 
original magnitude mjn. The efficiency matrix accounts 
for the incompleteness effects which are intrinsic to the 
detection algorithm, the flat-fielding and sky-subtraction 
errors. Depending on the way this matrix is computed, 
it also can include magnitude errors caused by crowding. 

The functional efficiency, E{mout), is usually defined 
as the fraction of sources detected with magnitude mout 
irrespective of their input magnitude, from the total 
number of objets that originally had magnitude mout^ 



E{mQut) 



A^det(mout) 
A'"orig(mout) 



(5) 



where A^dct(mout) is the number of sources which are 
detected with magnitud mout, while iVorig(mout) is going 
to represent the number of sources which originally had 
magnitude mout- Therefore, the functional efficiency is 
defined as a " detection rate" , and it is going to depend 
strongly on the initial number of input sources at each 
magnitude bin. 

From the previous definitions, the relation between 
E{mout) and P(min;mout) is given by 



and the functional efficiency, -©(mout), in the study with 
science frames; while we have only used the functional 
efficiency method in the study with synthetic images, 
because this is only used for a first estimation of the 
DETECT_THRESH range. As it will be commented later, 
the matricial method became very unstable at faint mag- 
nitudes when used to correct star number counts for ef- 
ficiency (see H5.4|l . Thus, we decided to apply the defini- 
tive efficiency corrections using the functional efficiency 
E{mout) computed with the study with science frames. 
Nevertheless, the efficiency matrices were essential when 
we performed our variation of the reliability analysis by 
CH03, as described in ^FT^ 

Once we computed P(min;mout) and i?(mout) using 
the science frames, the efficiency corrections were ap- 
plied as follows. We obtained the initial catalogues by 
running SExtractor onto the final U and B images. No- 
tice that the efficiency corrections must be applied to the 
number counts obtained from catalogues after removing 
spurious sources. This is because the efficiency matrices 
and functions are computed considering only the number 
of sources which are recovered from an original input dis- 
tribution, without counting the number of spurious detec- 
tions. Then, after rejecting the spurious sources in both 
catalogues as described in t|5.2.2l we proceeded to count 
sources by magnitude intervals according to the 3 size 
groups defined above. Differential number counts can be 
corrected for completeness by the following two methods: 

1. Using the efficiency function i?(mout)- Efficiency- 
corrected number counts per magnitude bin and 
unit area are given by 



5 (mo 



0.5-A-Es{mont)' 



(7) 



where we have defined magnitude bins of 0.5 mag; 
s = 1,2,3 refers to the size group; -/Vorig,s(wout) 
are the efficiency-corrected number counts for size 
group s and magnitude bin mout; Adct,s(wout) are 
the detected counts for that size group and mag- 
nitude bin corrected for spurious detections; A is 
the area over which we have made galaxy number 
counting; and -E's(mout) is the functional efficiency 
for the same size group and magnitude bin. Finally, 
contributions of each size group are added up to 
give completeness-corrected total number counts: 



£'(mout) 



J^Ym- ^dct(min,mout) Z^Vm ^(™in ; mout ) iVorig (min) 



A"orig(mout) 



A"orig(mout) 

(6) 

Note that if and only if A'orig(min) is constant for all min, 
then Eijnout) = Z]vmi„ -P("^in; mout)- 

The functional efficiency, _E(mout), is used more of- 
ten than the efficiency matrix P(min;mout) in number- 
count studies jl ladovich ct al. 2004; Metcalfe et al. 200^ 
iBershadv et alJll998j) . due to the instabihty which usu- 
ally arises when inverting P(mi„; mout) for applying the 
efficiency correcti on. Nevertheless, as pointed out by 
iHogg et al] l)1997|) . the advantage of using P{mi^; mout) 
is that it corrects not only for completeness errors due to 
the loss of sources, but also for photometric errors. 

We have corrected number counts for efficiency using 
the two methods: the efficiency matrix, -P(min; mout), 



A^orig .total (mout) 



A^o 



5 (mo 



(8) 



2. Using the efficiency matrix P{min', mout)- From the 
definition of -P(min; mout), the number of sources 
detected with magnitude mout irrespective of their 
input magnitude without having spurious detec- 
tions into account, is given by 



A"dot,s(mout (min) • ^'s (min; mout), (9) 

Vmin 

where A"orig,s(min) represents the source num- 
ber originally injected at m^n for the size group 
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s; -Ps ( Win ; TOout) is the efficiency matrix element 
P(TOin; TOout) for size group s; and Ndct,s{mout) are 
the detected counts for that size group and magni- 
tude bin corrected for spurious detections, as above. 
Thus, the efficiency-corrected number counts per 
magnitude and unit area are 



Norig.simnY = 0~5~A' I ^ ^dct,s(™out) ' Pg \wout,"lin) 
I Vrriout 

(10) 

We have defined (?7iout, "lin) as the element 
(wouti "T-in) of the inverted matrix of P{mi^; ruout), 
for objects from the size group s. Finally, contri- 
butions from the three size groups are added up to 
obtain total efficiency-corrected number counts at 
each band, using equation 

Error estimation for both methods is described in de- 
tail in Appendix ^ Briefly, when using the efficiency 
function, the _B(TOout) errors are quadratically added to 
those from counting statistics, while errors using the effi- 
ciency matrices P(min; TOout) are quantified by estimat- 
ing how number counts would change if maximum errors 
from counting statistics and from efficiency matrices are 
separately considered in equation (|10|) . 

5.1. Efficiency and Reliability Analysis for Synthetic 

Images 

As commented in the previous section, the study of 
the efficiency and reliability using synthetic images must 
be interpreted with care, due to the difficulty in repro- 
ducing the sky conditions of the image and the profiles 
of simulated sources, which usually are regular in excess. 
So, we must remark that we have used this analysis us- 
ing synthetic images only to narrow down the range of 
DETECT_THRESH values for the definitive study with the 
science images, which is more detailed and trustworthy. 

Using the art data package in IRAF, we created these 
artificial images with Poisson background noise of same 
RMS as our real images. Synthetic disk galaxies and 
stars of known magnitudes and sizes were added at given 
positions. Magnitudes span the range for stars and galax- 
ies in our science frames. The number of input galax- 
ies at each magnitude bin was chosen to reproduce an 
initial guess for the slope of the counts in each band, 
while the number of stars was the estimated using the 
iBahcall IT Soncira (198l3 model for the Galactic coor- 
dinates of the GWS. The sizes of all the stars and of half 
of the galaxies were set to the PSF FWHM in each band, 
while the other half of the synthetic galaxies was twice 
bigger. Detection efficiencies and spurious fractions as a 
function of source magnitude were determined by run- 
ning SExtractor with different DETECT_THRESH values 
(sec Radovich ct al. 2004; Metcalfe ct al. 2001; Lin ct al. 
119981 for more details of this procedure). In Figure 0] 
we show the efficiency functions, £^(mout)7 and spurious 
fractions for several DETECT_THRESH values, computed us- 
ing the artificial B image of the GWS, as defined in 
SjSl The magnitude of 50% efficiency (to5o%) increases 
as we lower the detection threshold, while spurious de- 
tections quickly rise at magnitudes approaching TO5o%- 
From Figure 0] and the corresponding distributions for 



the U band which behave similarly, we concluded that 
detection thresholds ranging from 0.4 to 0.7 provided a 
good compromise between high detection efficiency and 
low spurious fraction in both bands. Thus, we decided to 
use the following DETECT_THRESH values for the definitive 
study of incompleteness using the science frames: 0.4, 
0.5, 0.6, and 0.7. 

Notice that spurious fractions can be easily determined 
I when using synthetic frames, because all the sources in 
rthe artificial images have been inserted by us. But this is 
' not the same when using science frames, where we have 
a mixture of those sources we insert and of those which 
were there originally. In fact, the spurious analysis when 
using science frames is more complex (see t|5.2.2|l . 

5.2. Efficiency and Reliability Analysis for Real Images 
5.2.1. Efficiency 

We have carried out an extensive series of simu- 
lations to quant ify the detection efficiency using the 
real images (see iHog g et all 119971 iHuane et al. 2001^ 
IKiimmel fc WagneJ 2001. CH03, among others, for a de- 
tailed description of the procedure) . Firstly, we found the 
brightest source present in our images, in each one of the 
three size groups defined in ^ Then, these three selected 
sources were inserted several times in the science frames 
following a flat magnitude distribution (18< B (or U) 
<28) at random locations. No constraints were imposed 
on the source positions, in order to include the effects of 
source confusion in the computation of P(min; mout) and 
E{mout)- This process was repeated 50 times for each 
size group, injecting 2,000 sources in each simulation, for 
a total of 300,000 objects. SExtractor was run on each 
simulated image with DETECT_THRESH=0.4, 0.5, 0.6, and 
0.7. For each simulation, we computed the efficiency ma- 
trix, P(min', ?Tiout), and the functional efficiency, E{mout) 
(see their definitions in SjSJ. 

In Figure El we show the average i?(r7iout) from all the 
simulations, as a function of magnitude and size group, 
for the four selected values of DETECT_THRESH, in both 
B and U bands. Error bars ar e the quadratic addition 
of the RMS of all simulations ijBevingtonI 1 1 969(1 . In all 
the cases, E{mout) shows a gentle decline, followed by 
an abrupt drop near the detection limit. The 50% ef- 
ficiency magnitudes for a given object size differ by as 
much as ~0.6 mag when changing DETECT_THRESH; while 
they are ~0.7 mag deeper for point sources than for large 
objects of the same magnitu de. These trends are si milar 
to those found by CH03 and lBershadv et all l|1998|) . We 
have also found that our synthetic frames overestimate 
the efficiency of SExtractor on to t he science ima ges, a 
fact that corroborates what iBershadv et al.l ((,1998) and 
CII03 reported. It can be noticed just comparing the 50% 
efficiency magnitudes obtained using synthetic frames for 
DETECT_THRESH=0.5,0.7 in B (see the Figure|l|), with the 
achieved ones using the science frames (see the Figure^}. 

The sizes of the selected objects were approximately 
typical for their size groups, except for the object rep- 
resenting the "intermediate-sized" objects in B, which 
is large into the range of its size bin. Nevertheless, we 
have checked that this fact does not underestimate the 
typical detection efficiency of its size group for magni- 
tudes less than the 70% efficiency magnitude for a fixed 
DETECT_THRESH. For higher magnitudes, the curve of the 
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Fig. 4. — Comparison of detection efficiencies i?(m,out) and spurious fractions as a function of the detection magnitude rriout, obtained 
using the synthetic B image for several DETECT_THRESH values. Drawn curves include all simulated objects, irrespective of their size. The 
situation at the U synthetic image is completely analogous. Upper panel: Efficiencies as a function of the magnitude of detection. Lower 
panel: Fractions of spurious sources detected at each magnitude interval relative to the total number of sources inserted at this magnitude. 



efficiency function would be displaced < 0.15 mag to 
fainter magnitudes. But the changes in the results would 
be negligible in the last significant bin of magnitude (the 
50% efficiency magnitude bin), because we have defined 
bins of 0.5 mag and the effciency drops from 70% to 
50% in less than 0.5 mags (see Figure EJ. Therefore, we 
can consider that the surface brightness of each selected 
brightest object represents the typical surface brightness 
in its correponding size group. 

B-band average P(mi^^; rriout) values for point-like ob- 
jects are plotted in FigurelHlfor the four detection thresh- 
old values we studied. Errors associated to these matrices 
are also shown in the figure; they have been computed in 
the same way as the i?(r7iout) errors. Similar behaviours 



are seen for the three size groups, and for both filters. 

As we have commented above, we corrected num- 
ber counts for incompleteness effects using the effi- 
ciency matrix, P(min; mout), and the functional effi- 
ciency, i?(mout)- The former became very unstable 
at faint magnitudes when used to correct star number 
counts for efficiency, as described in ii5.4l So, we decided 
to apply the definitive efficiency corrections using the 
functional efficiency E (rriout) only. The efficiency matri- 
ces were useful when we performed our variation of the 
reliability analysis by CH03 (see H5.2.2|l . 

We must remark that the P(min; ?nout) determination 
does not depend on the used input source distribution, 
and thus, it can be used for correcting for completeness 
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Fig. 5. — Left panels: Detection efficiency functions, -E(mout), given as the ratio of detected to input sources as a function of detection 
magnitude mout, for the three apparent size groups in the B image, using DETECT_THRESH=0.4, 0.5, 0.6, and 0.7. Errors are plotted for 
DETECT_THRESH=0.6 only, because they are similar for all the detection thresholds. Each line style corresponds to one of the DETECT.THRESH 
values according to the legend. The 50% efficiency magnitudes for each size group are listed in each frame. Right panels: The same for the 
U image. 



whatever the real input distribution is. But this is not 
true for i?(mout), which would have need an input source 
distribution as similar as possible to the real one. Nev- 
ertheless, this is not feasible when using science frames, 
because these images are already very populated by real 
objects. An aditional population of inserted objects as 
large as the original one would crowd the image exces- 
sively, and as a result, computed efficiency matrices and 
functions would overestimate the effects of source con- 
fusion. Thus, an additional error arises when number 
counts are corrected using an efficiency function which 
has been computed with an input source distribution flat 



in magnitudes (instead of an input one that reproduces 
the initial slope of the counts). The efficiency function 
error that arises at each magnitude bin can be estimated 

as _ 

AE{mout) = 77- ^ , (11) 

where i?(mout) is our efficiency function computed using 

the flat input distribution of sources, and i?(mout) rep- 
resents the efficiency function that would be obtained 
using an input distribution with the initial slope of 
the counts (i.e., the correct one). We have estimated 
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JTiout) of finding an injected source with 

initial magnitude at output magnitude mout, for point-like objects in the B image, using DETECT_THRESH=0.4, 0.5, 0.6, and 0.7. Solid 
line indicates mBi^ = mBout- Right panels: Standard deviation matrices associated to the efficiency matrices at the left. Similar patterns 
are found for all size groups and for the U image. 



E{mout) at each magnitude bin as foUows: let us con- 
sider an input source distribution identical to the one we 
detect once the spurious sources have been subtracted 
for each size group s, A'^siopc^sCmout)- As P(min;mout) 
is independent of the used input source distribution, 
the correct distribution of sources we are going to de- 
tect is A'"dot,s("T.out) = Z^Vmin ^slopo,s("lin)-Ps(min; mout)- 



Moreover, if we had used the correct efficiency function 
E{mout), the detected source distribution would have 
been the same: A'^det^(mout) = ^siopo,s(mout) ■ ^^(TOout)- 
So, we can deduce Eimout) from the two previous ex- 



12 



Eliche-Moral et al. 



pressions: 

~. , Evmi„^slopc,4win)Ps(TOi„;mout) 

E{mout) = 7- — T , (12) 

We have estimated the error that arises in using E{mout) 
instead of £'(mout) through equations p2(l and (|1H) . 
and considering the detected source distribution once 
the spurious have been removed as the input distri- 
bution A'siopc,s(TOout)- The maximum error Ai?(mout) 
reached in both bands (and thus in the number counts) 
is < 10% for magnitudes m < 77159% i^i each size group, 
which is less than the total error from statistical counting 
and efficiency errors using the computed E{mout)- Nev- 
ertheless, the E{mout) determination has itself a high 
uncertainty which would overestimate efficiency errors, 
and that arises from the division by A'siopc,s(™out) in 
those magnitude bins where we have low statistical sig- 
nificance, and from the accumulation of the errors from 
Psimin, rriout) and A^siopc,s(mout) due to the error propa- 
gation of equation H12() . Therefore, we decided to ignore 

this error and use E{mout) instead E{'mout) for correct- 
ing number counts. 

5.2.2. Reliability 

Reliability has been characterized us ing the method 
descri bed in CH03, based on that used bv lBershadv et al.1 
()1998ft . It consists on creating two half-time exposured 
images from exclusive halves of the data. As spurious 
sources are basically noise correlated structures or peaks 
that appear in one of the subexposures, they must appear 
only in one of the two half-time images, the one that 
was built including the correponding subexposure. By 
running SExtractor in double-image mode, photometry 
on both half-time images is measured at the positions of 
sources detected by SExtractor in the total-time image. 
Those spurious sources detected in the total-time image 
will have very different magnitudes when measured in 
the half-time images, as they have a high probability of 
appearing only in one of the two half-time images; while 
the real sources will have very similar magnitudes when 
measured in the two half-time images, because they have 
a high probability of appearing in both half-time images 
with approximately the same flux. We must notice here 
that this method could not be applied to B data directly 
due to their dithering sequence (see Table 01 . However, 
the results obtained through this method in the U band 
can be used for characterizing reliability in the B image, 
as we will comment below. So the {S/N)iiyn method for 
characterising SExtractor reliability was carried out on 
the U image only, and extrapolated to the B image at 
the end. 

In order to identify spurious detections, CH03 used a 
signal-to-noise ratio {S/N) criterion. Hereandafter, we 
will define the S/N of a source as the ratio between its 
flux and its flux error measured by SExtractor. The 
method consists on assuming that all the sources with 
S/N below a given limit, (S'/A^)ii,n, in one or in both 
of the half-time images are considered as "false detec- 
tions", and hence, rejected from the catalogue. Never- 
theless, into this group of sources labelled as "false detec- 
tions" , we will be rejecting sources that actually are spu- 
rious detections (henceforth called "rejected truly spu- 
rious"), and sources that in fact are real objects whose 



S/N's also obey the imposed criterion (called "rejected 
real sources", hereandafter). As we increase the value of 
the imposed (S'/iV)iini, the probability of being rejecting 
more truly spurious sources increases. But, at the same 
time, the probability of being rejecting a real source as 
a "false detection" also raises. There will be a limiting 
{S/N)[^^ equal to the one of the highest S/N that the 
truly spurious sources exhibit in the catalogue, so that it 
rejects most of them from it. Using {S/N)ii^ > {S/N)[^^-^ 
will not reject more truly spurious, because the major- 
ity of them have already been rejected with the lower 
{S/N)[^^, but it will reject more real sources as "false 
detections" (those with {S/N){.^ < S/N < iS/N)ii^). 
This is due to the fact that the population of real sources 
has higher S/N on average than the truly spurious pop- 
ulation, which is mostly related to the sky noise. There- 
fore, the fraction of rejected truly spurious sources will 
become constant for values (S'/iV)iii„ > {S/N)[^^ in each 
DETECT_THRESH, and this constant fraction can be consid- 
ered as the maximum fraction of truly spurious sources 
in the catalogue. 

For summarizing, for a fixed DETECT_THRESH, low val- 
ues of the {S/N)iim will reject a low number of sources 
as "false detections" from the catalogue, but with a high 
confidence level in that the majority of them will be truly 
spurious instead of real sources. On the other hand, 
high values of the {S/N)ii^ will reject a high number 
of sources as " false detections" , but we will have an un- 
known mixture of real and truly spurious sources into the 
rejected group. For {S/N)rnn values higher than an un- 
known limiting {S/N)[^^^, the majority of the truly spu- 
rious will be labelled as "false detections" (so the num- 
ber of truly spurious which are rejected will be approxi- 
mately constant), but the number of real sources labelled 
as "false detections" will keep increasing. 

Moreover, when a good discrimination between real 
and truly spurious sources is achieved for a given 
{S/N)iim, the majority of sources which are rejected as 
"false detections" will be truly spurious. Therefore, they 
will have very different magnitudes when measured in the 
half-time images, because a false detection has a negli- 
gible probability of being detected in both half-time im- 
ages with similar magnitudes. If we plot the histogram of 
Am = mhaif.i — 7nhaif,2 for the rejected sources with this 
fixed (S'/iV)iim, it will show a bimodal structure, as the 
majority of the rejected sources are truly spurious. On 
the other hand, if the fixed (S'/A^)iim does not distinguish 
properly between real and truly spurious sources, there 
will be two populations mixed in the Am histogram: the 
rejected truly spurious sources, which will draw the bi- 
modal structure, and the rejected real sources, which will 
contribute to a narrow single-peaked distribution cen- 
tered in Am — 0, because the real sources have a high 
probability of having very similar magnitudes in the half- 
time images. 

In Figure [7| the histograms of magnitude differences 
between the half-time images for candidate "false de- 
tections" in the U band are drawn with solid lines, for 
different (S'/iV)iini and for the four DETECT_THRESH val- 
ues we decided to study (see M5.1(l . Seventeen values 
ranging from (S'/iV)iii„=0.2 to (S'/iV)iii„=10.0 were stud- 
ied, but we have plotted results which correspond to 
(5'/A^)iim=1.4, 2.2, 2.8, 3.5, and 4.5 for clarity. When 
a good discrimination between real and spurious sources 
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Fig. 7. — Histograms of magnitude differences between the half-time images for candidate "false detections", and their corresponding 
estimated real source distributions. Solid lines: Histograms of magnitude differences between the half-time images for candidate "false 
detections" with different {S'/A'')iini and for the DETECT.THRESH values we are analizing. Corresponding (S/A'')iini are 1.4, 2.2, 2.8, 3.5, and 
4.5, according to the legend in the figure. Dashed lines: Estimated real source distributions compatible with each one of the the rejected 
source histograms. 



is achieved, the histograms show a bimodal structure, 
as for the lower values of (S'/A^)iini. Conversely, as we 
raise (S'/A^)iim, sources with higher S/N, and hence with 
more probability of being real, are classified as "false 
detections". This situation is responsible for the single- 
peaked, centered distributions of Am = in the cases 
of high (S'/A^)iim. As can be seen from the figure, select- 
ing the combination DETECT_THRESH-(S'/iV)ii„i from these 
histograms is tricky, because histograms evolve from dou- 
ble to single-peaked gradually; i.e., there is always a pop- 



ulation of rejected real sources mixed up with the distri- 
bution of rejected truly spurious detections, considered 
as "false detections" by this method. Until now, the se- 
lection DETECT_THRESH-(S'/iV)iim was subjective, depend- 
ing on which histogram seemed to be more double-peaked 
and populated at the same time. 

We have developed a self-consistent procedure for 
choosing the most adequate DETECT_THRESH-(S'/A^)iim 
combination in an objective way, in order to minimize 
the number of real sources which are rejected as "false 
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detections" , and to maximize the number of rejected 
truly spurious sources and also the number of real de- 
tections that remain in the catalogues. The proce- 
dure is described in detail in [Elichc-Moral & BalccUs 
l|2005. hereafter EB05). so we briefly summarize neces- 
sary information for this paper. The method consists 
on estimating statistically the population of real sources 
which is mixed up with the one of truly spurious sources 
in the histograms of "false detections" of Figure [7| For 
computing it, the following justified assumptions are used 
in the U band: 

1. The detected source distribution directly obtained 
by SExtractor, N^cti'm): was taken as a first ap- 
proximation of the distribution of detected real 
sources, iVdet,reai(w), being m the detection mag- 
nitude. 

2. The distribution iVdet,reai('7i) was truncated for 
magnitudes mil > 26 mag, because 50% efficiency 
magnitudes are less than mil = 25 mag for all the 
DETECT_THRESH values we are analysing. Therefore, 
a source with mU > 26 will have a ~ 100% proba- 
bility of being a truly spurious detection, and these 
magnitude bins will not contribute to the real pop- 
ulation in the histograms of magnitude differences 
of Figure 13 

3. We have also considered that all the sources at 
the histograms of AmU that exhibit AmU — 
wJ7out,i — wJ7out,2 = are real detections. This 
is justified by the fact that the Poissonian proba- 
bility of detecting two spurious sources at the same 
position with the same magnitude at both half-time 
images is negligible. 

Once the distribution of real detections compatible 
with each histogram of magnitude differences is com- 
puted for each (S'/iV)ii,„-DETECT_THRESH combination, we 
can easily estimate the fractions of real sources and truly 
spurious we are rejecting from the total number of de- 
tections, as well as the fraction of truly spurious sources 
that remain in the catalogue in each case. Notice that 
we can estimate also the limiting {S/N)[i^ so that using 
{S/N)\i^ > {S/N)[i^ will not reject more truly spurious. 
This will be the {S/N)ii^ value from which the number 
of rejected truly spurious remains constant although wc 
raise it. 

In Figure 13 we have overplotted with dashed lines the 
estimated distributions of real sources compatible with 
each histogram of magnitude differences in the half-time 
images. As (S'/A^)iim increases, the estimated real distri- 
bution approaches the corresponding histogram, because 
most of the rejected "false detections" are real sources. It 
is remarkable that the widths of these real source distri- 
butions (FWIIM-0.4-0.7 mag for aU the DETECT.THRESH- 
(S'/A^)iim combinations) are consistent with the photo- 
metric error for sources at those S/N on the half-time 
images, a fact that supports their robustness. 

In Table im we have compared the estimated fractions 
of rejected real sources, rejected truly spurious sources, 
and truly spurious that are non- rejected (i.e., that remain 
in the catalogue), for the four values of DETECT_THRESH 
and five of the seventeen values of {S/N)\ini wc arc 
analysing in the U band: S/N=1.8, 2.2, 2.5, 7, and 10. 



The procedure for getting the values from Table the 
arguments to set the (S'/iV)iii„-DETECT_THRESH combina- 
tion of values, and the results are given in detail in EB05. 
Thus, we will just explain them briefly here. 

Firstly, notice that the fraction of truly spurious 
sources becomes constant for {S/N)iiyn = 6, 10 in 
each DETECT_THRESH in the Table. In fact, it is con- 
stant for {S/N){.^ > 2.2 in DETECT_THRESH=0.4, for 
[S/NYy^^^ > 3.0 in DETECT_THRESH=0.5, for iS/N){.^^^ > 
3.5 in DETECT_THRESH=0.6, and for {S/N)[.^^^ > 4.0 
in DETECT_THRESH==0.7. This means that the maxi- 
mum number of truly spurious that are rejected in each 
DETECT.THRESH is reached with the corresponding S/N(.^ 
value. Then, no advantage arises in using (S'/A^)iim > 
{S/N)[^^^. We have also found that the maximum 
fraction of truly spurious that are rejected for all the 
{S/N)iim is the same for DETECT_THRESH<0.6 29%), 
but different for DETECT_THRESH=0.7 25%), as you 
can infer from Table ITTl Therefore, DETECT_THRESH=0.6 
seemed to be the limit between two different behaviours 
in detection. Although the total fraction of rejected 
truly spurious (respect to the total number of detec- 
tions) is higher if we use DETECT_THRESH<0.6 than 
using DETECT_THRESH=0.7, it can be controlled, and 
it is possible to reach deeper magnitudes than using 
DETECT_THRESH=0.7. With an adequate choice of the 
{S/N)ii^ (see Table [TTll . we can get similar fractions of 
truly spurious rejections, of real source rejections, and 
of truly spurious remaining in the catalogues, and in- 
crease the number of detections. Thus, we discarded 
the value DETECT_THRESH=0.7. Moreover, {S/N)ii^^2.2 
gave a good compromise between rejected real and false 
detections for DETECT_THRESH<0.6. So we finally de- 
cided to use DETECT_THRESH=0.6 and (5/7V)iii„=2.2 for 
the U filter, and hence, the 50% efficiency magnitude for 
point-like sources is J7DET=o.6(Eff=50%)=24.83 mag (see 
Figure El . This selection allowed us to reject ~22% of 
truly spurious sources in the catalogue (see the values 
in Table [TT|l . As the maximum fraction of truly spuri- 
ous sources for DETECT_THRESH=0.6 was estimated to be 
~29% of the total number of detections, then its is easily 
deduced that approximately ^7% of the final catalogue 
are truly spurious detections, the bulk of them at mag- 
nitudes fainter than the 50% efficiency magnitude. With 
this choice, ~3% of the real sources are rejected from 
the catalogue. Using (S'/iV)ii,„=2.2 in half-time images 
for rejection, sources considered as real detections have 
roughly {S/N)\i^ ^3.1 in the total-time image. 

As commented previously, this method could not be 
applied to B data directly because it was not possible to 
obtain two half-time images from the dithering sequence 
(see Table 0J). If we would have created two complemen- 
tary images without using all the subexposures in the 
B band, these two images would have been exposed less 
than itotai/2. Therefore, we would have had a lower con- 
fidence in detecting spurious sources than in the U band, 
and we also would have lost the spurious sources corre- 
sponding to those subexposures not used for creating the 
two complementary images. Nevertheless, the results for 
the U fllter can be used for estimating the best combi- 
nation of DETECT_THRESH-(S'/iV)iii„ for the B image as 
follows. The B and U images had similar flux charac- 
teristics, source distributions and Poissonian sky noises, 
and hence, the detection threshold for U image should be 



U and B Number Counts in the Groth Strip 



15 



valid for the B image too. This is because, although a^ky 
is higher in the B image, B is deeper than U, and their 
efficiency matrices and functions behave similarly (see 
Figures [lland|n|). Therefore, we use DETECT_THRESH=0.6 
for the B-band also, and obtain a 50% efficiency magni- 
tude of BDET=o.6(Eff=50%)==25.46 mag (see Figure EJ. 
In order to establish the (S'/A^)iim for the B image, we 
have considered that detections only depend on the S/N 
of the sources once the minimum area has been fixed. 
When the S/N of a source is below the fixed (S'/A^)iim, 
its flux per unit time and area must be lower than the 
flux that corresponds to that (S'/iV)iim. In EB05, the 
equations that relate the flux of a source and its S/N are 
shown as a function of several image parameters. This 
flux can be expressed in units of the sky sigma of the 
image. Therefore, a source with (S'/A^)iim=2.2 in the U- 
band with the minimum area and for t — Uotai/'^ would 
have a flux per pixel and per unit time equal to 1.73 times 
the sky sigma of the total-time U image. Extrapolation 
of this result to the B image is justifled due to the simi- 
larity of U and B images, and hence we have considered 
that a source in the _B-band will only be considered a real 
detection if its flux per pixel is greater than 1.73 a sky of 
the B image, which implies a (S'/A^)iim=2.8 for a source 
with the minimum area in the B total-time image. 

Finally, we used the identified "false detections" in the 
U band with the selected combination of DETECT_THRESH- 
{S/N)ii^ for correcting the corresponding catalogue (ex- 
tracted using the selected DETECT_THRESH=0.6) for spu- 
rious detections. Hereandafter, we are going to call 
"spurious sources or spurious detections" to all the 
"false detections" that have been rejected. As the same 
method can not be applied to the B band, we cor- 
rected for spurious detections the B catalogue obtained 
with DETECT_THRESH=0.6 rejecting all the sources with 
(S'/A^)iim < 2.8 in the B total-time image. 

Our procedure for estimating the distribution of real 
sources in the Am histograms is able to fix objectively 
the DETECT_THRESH-(S'/7V)iim combination which mini- 
mizes the number of rejected real sources, maximizing 
the number of true spurious rejections, and reaching the 
deepest limiting magnitude at the same time. In fact, 
the subjetive selection of {S/N)\iyn by CH03 leads them 
to remove the 13.7% of the sources with S/N > 5 in 
their data, while we have rejected < 1% of sources with 
S/N > 5, < 9% of sources with 3 < S/N < 5, and - 13% 
of sources with S/N < 3. Our method allows us to esti- 
mate that only ^--^3% of real sources with S/N < 3.1 are 
rejected from our catalogues. 

5.3. Galactic Extinction 

Even at the high Galactic latitude of the GWS field 
(6=60°), Galactic extinction affects U- and i?-band num- 
ber counts. Lack of extinction corrections probably ex- 
plain the differ ences between published £/-band num - 
ber count data l|Heidt et a,].ll20Mm,a,dnvich et a,].ll200l . 
We have computed Galactic extinction corrections for 
each source, i n both the U and B filters, using the 
ISchlegel et all l|199at) extinction maps . The average 
value is E{B -V) = 0.011037 ± 0.000004, ranging from 

Dust maps and software for computing the extinction correc- 
tions for eac h source were downloaded from : 
http: / / astron.berkeley.edu / davis / dust /local /local, html | 
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Fig. 8. — Star counts in GWS field derived using the DEEP cata- 
logue, corrected for detection efficiency using "functional" method 
(squares) and "matricial" method (triangles) for both bands. Di- 
agonal crosses are raw star counts in each filter. Solid line is the 
star counts prediction at each band for the GWS {b = 60°, I = 95°) 
from the Bahcall & Soneira (1980, 1981a) Galaxy model (see text 
for more details). Error bars correspond to counting statistics, 
added in quadrature to uncertainties from the efficiency correction. 



0.008 < E{B -V)< 0.018, which translates into a max- 
imum Galactic extinction correction of An 0.11 mag 
and Ab ~ 0.08 mag. 

5.4. Star-galaxy separation 

Number counts for "point-like" sources (defined in 
^5.2|) were corrected for star counts in order to obtain 
galaxy number counts. Star identifications were directly 
obtained over the area of overlap between our images 
and the HST Groth survey by cross-correlating our cat- 
alogue with the F606W Medium Deep Survey (MDS) 
catalog of the GWS field^s, using the MDS star iden- 
tifications (Ratnatunga et al. 1999). The resulting star 
counts were corrected for detection efficiency using the 
efficiency functions and matrices (see i l5.2.1|l . scaled to 
the total areas for U and B, and subtracted from the 
efficiency-corrected "point-like" number counts. Star 
number counts in U and B are shown in Figure |H1 
corrected for efficiency errors using the matricial and 
functional methods. Both procedures give similar re- 
sults for our entire magnitude range, but the inversion 
of the efficiency matrices became strongly unstable for 
B{U)>24: mag. Therefore, we decided to use the func- 
tional efficiency method for correcting final galaxy num- 
ber counts. We have also compared our star count mea- 

1^ M edium Deep Survey, MDS: 
|http: //archive .stsci.cdu/mds 
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TABLE 5 

U Differential Star Counts in the GWS Field 



u 


^ star, raw 


Wstar.l 




ci.l 






0"I,2 


(1) 


(2) 
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(8) 


18 25 
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18 75 






900 90 


139.78 
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212.70 


153.56 


19 25 
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246.81 


134.46 


65.15 
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314.90 
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239.36 


181.48 


20.25 
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zvj. 1 O 
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168.52 




9*37 Qft 
ZO 1 .ou 


1 77 


21.25 


404.64 
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205.63 


143.64 


365.02 


198.07 


135.82 


21.75 


657.54 


682.67 


247.63 


187.18 


725.96 


253.45 


193.00 


22.25 


708.12 


741.54 


256.88 


196.22 


595.24 


239.57 


179.28 


22.75 


1315.08 


1395.82 


332.97 


273.78 


1588.32 


362.57 


300.18 


23.25 


1719.72 


1842.21 


376.94 


318.20 


1257.36 


282.47 


235.81 


23.75 


2174.93 


2395.14 


432.02 


372.85 


4076.67 


694.77 


600.13 


24.25 


2377.25 


2713.88 


469.13 


408.58 








24.75 


2225.51 


2927.68 


535.40 


466.94 








25.25 


1466.82 


3990.05 


1011.32 


881.00 








Note. — 


Col. (2) are U raw 


star counts 


in units of N mag 




(once the 


spurious 



where removed from the catalogue); Col. (3) shows efficiency-corrected U star counts by the 
functional method in units of N mag"""" deg~^ (denoted by subindex 1); Col. (4) and (5) are 
the 1-cr confidence upper and lower errors associated to the functional-efficiency method; 
Col. (6) shows efficiency-corrected U star counts using the matricial method in units of N 
mag"""" deg~^ {denoted by subindex 2); Col. (7) and (8) are the l-cr confidence upper and 
lower errors associated to the matricial-cfficiency method. 



TABLE 6 

B Differential Star Counts in the GWS Field 



B 
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.^star.raw 
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Note. — 


Columns are as in 


Table El 

















sureme nts to the predictions from the lBahcall &: Soneiral 
((1981a) model of the Galaxy for the coordinates of the 
GWSi^ (see alsolBahcall & Soneira , 1980i Bahcall 1986; 
iRatnatunga & Bahcall 198 51)- Adopted Gal axy parame- 
ters are from T able 1 o fiCaban ac et alJ lll998D , LF param- 
eters are from lMamon fc Soneiranil982tl. and r epresen - 
tative color terms are from iBahcall fc Soneiral l|1981bj) . 
Predicted counts from the model are overplotted in Fig- 
ure |H1 Model and measurements agree in both bands at 
intermediate magnitudes. Divergen ces at the f a int end 
are to be expected, as the Bahcal Pfc Soneir al i 1981 all 
model is only reliable down to £? = 20 ( Bahcall^^e^^ 
[1985: Laske r et al . 1987; Santiago et al. 1996; Bat h et al. 
11996ft . At the bright end, the models predict ~10 sources 

The source code for computing the star counts for different 
filters and Galactic coordinates is kindly available by the author in 
the Astrophysics Source Code Library Archive, BSGMODEL: The 
Bahcall-Soneira Gala xy ModeL 
jhttp : / / ascl . net/bsgmodel . html j 



per 0.5-magnitude bin at f/(or B)^20 in the area of the 
HST GWS data. The divergences here are likely due to 
low number statistics in our star count measurements, 
probably due to the rejection of saturated objects in the 
catalogues. For the bright regime between C/(or i3)~20 
mag and the saturation limit^^, we have resorted to 
SExtractor's stellarity parameter (CLASS.STAR > 0.9 
for both U and B), shown by Capaccioli et al. (2001) to 
reliably identify stars at bright magnitudes in blue filters. 

Results from star counts in U and B are listed in Ta- 
bles |S1 and El Columns indicate raw star counts with 
the spurious sources subtracted; and efficiency-corrected 
star counts using functional efficiencies and matrices, to- 
gether with the estimated upper and lower limits for both 
methods. The correction in log(iV) due to stars is lower 
than 0.05 index at fainter magnitudes. Errors in star 

1'^ A point source with FWHM=1.3"saturatcs at B ~20.0 and 
U ~18.3 mag in our images 
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Fig. 9. — Galaxy number counts from GWS data in IJ and 
B over ~900 arcmin^. The comparison wit h previous works in- 
cludes only tabulated, CCD-based counts (see'Gardner et al.'igQfT; 
l^lonteri ct al. 2000; Metcalfe ct al. 2001; Yasuda et al. 2001; 
IHuang et al...2001b.: .Kiimmel fc W agner 2001; Capak et al. 2004). 
The faintest point in each band corresponds to the 50% efficiency 
magnitude, while the brightest end reaches the saturation limit of 
the INT/WFC in each band. 



counting are computed in the same way as the errors 
in galaxy counting. Lower and upper errors from star- 
counting will be cuadratically added to the correspond- 
ing upper and lower errors that arise from efficiency cor- 
rections and statistical counting to obtain final number 
counts errors (see Appendix IX)) . 

6. RESULTS AND DISCUSSION 

6.1. Differential number counts 

The final source catalogues were obtained by run- 
ning SExtractor on the U and B images with 
DETECT_THRESH=0.6, excluding s purious candidates ac- 
cording to the criteria outlined in t|5.2.2l Extinction cor- 
rections were applied to all the sources in the catalogues 
(see H5.3|l . For the number count study, the edge areas 



with lower exposure due to the dithering pattern were 
excluded. Final useful areas for number counts were 888 
arcmin^ for B and 846 arcmin^ for U (see Table 0)). The 
division into the three size groups defined in ^jSjallowed us 
to correct for completeness in each size group using the 
functional efficiency method, as described in H5.2.f I For 
summarizing, final counts are obtained from raw counts 
through: Galactic extinction correction f i|E3)l : correction 
for spurious detections f H5.2.2|l : detection efficiency cor- 
rection using the functional efficiency method ft |5.2.1|l : 
and subtraction of the star counts from Tables|5l and |B1 

Our results for U and B galaxy differential number 
counts are summarized in Tables IT3 and [T51 respectively. 
Magnitudes are in the Vega system. Raw counts in the U 
and B bands are listed in Col. (2)-(4) in these Tables, for 
each one of the three source size groups defined in H5.ll 
Spurious corrected number counts are shown in Col. (5)- 
(7). Applied efficiency correction factors for each size 
group appear in Col. (8)-(10). Col. (11) gives the final, 
differential galaxy number counts TV per magnitude and 
per deg^, and Col. (12) and (13) fist the upper and lower 
la errors for N, computed as explained in Appendix IXI 
For convenience, we list log(iV) in Col. (14). Counts are 
derived from 0.2350 deg^ and 0.2467 deg of sky in U and 
B, respectively. The range of our counts is 18 < J7 < 25 
and 19.5 B < 25.5. The bright limit is set by satura- 
tion on our CCD frames. The faint limit is set by our 
50% detection efficiencies, whic h are 1/50% = 24.83 and 
^50% = 25.46 for point sources f H5.2.1|l . Very few sources 
are detected in the large size group (see Col. (4) of Tables 
ll2land [TS|l . and these are only found in the faintest mag- 
nitude bins. Visual inspection shows that most of these 
sources are single- or multiple-peak structures embed- 
ded in a region of higher background, which SExtractor 
takes as single extended sources. While they may include 
real sources (note that these detections have survived the 
spurious source filter), it is very unlikely that sources in 
the large size bin trace single, faint, extended galaxies. 
For consistency, we apply efficiency corrections to them, 
and add them to our galaxy counts. Their contribution 
to the total counts is entirely negligible. 

We plot U and B number counts in Figure El together 
with literature data. We have only included number 
count data coming from CCD observations which were 
given in tables, and convert U and B data from different 
photometric systems to Vega magnitudes in the Landolt 
system. Our counts are in excellent agreement with the 
other studies. Note that literature data near the bright 
and faint ends of our counts come from independent stud- 
ies, and that our data bridge the gap between large area, 
shallow surveys and deep, pencil-beam surveys. The dis- 
persion among the different authors, which is greater in 
the J7-band, is probably due to different completeness 
and spurious rejection corrections, and to the absence of 
Galactic extinction correction for the majority of stud- 
ies. Clustering and cosmic variance can also contribute 
to the dispersion present in the U and B number counts 
from the different authors. We have estimated that the 
contribution due to clustering fluctuation is ~0. 1-0.4 the 
statistical counting error s in the range 20 < U(B) < 24 
mag f or both filters (see ■Tones ct afigQl': "Met calfe et alJ 
[T995HVoIo"nteri et al. 2000: Yasuda et al...200ll and ref- 
erences therein). 
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TABLE 7 

U AND B Galaxy Count Slopes from Several Surveys 



Reference 



dlog(N)/dm 



Magnitude Range 



C/-band 



Koo (1986) 
^^iTTiams et al. 



Hogg et al. fl997l 
^Pozzetti et al. (1998) 
Montana et al. (1999) 
Volonteri et al. (2000) 
Metcalfe et al. (2001) 
Radovich et al. (2004) 
Caoak et al. (2004^ 
Present work 



~0.4 20.0 < C/(Vega) < 22.0 

~0.05 26.0 < F300Vy(AB) < 28.0 

~0.40 23.0 < F300Ty(AB) < 26.0 

~0.467 20.3 < C/(AB) < 26.3 

~0.49 20.0 < C/(Vega) < 25.0 

~0.49 20.0 < C/(Vega) < 25.0 

0.47±0.05 22.25 < F300Ty(AB) < 27.25 

~0.4 18.0 < C/(Vega) < 25.0 

0.54±0.06 18.5 < C/(Vega) < 22.5 

0.526±0.017 20.0 < U(AB) < 24.5 

0.48±0.03 21.0 < C/(Vega) < 24.0 



B-band 



Tvso n (1988) 
Jone s et al. (1991) 
^Metcalfe et al. ( 1991) 
'Metcalfe ct al. ( 199^ 
Williams et al. (1996'l 



Bertin fc Dennefeld (1997) 
Pozzetti et al. (1998) 

^^^^^^^^^T ^^^^ 

Volonteri et al. (2000) 



Metcalfe et al. (2001) 
Huang et al. (2001^ 
'Kummel fc Wagner (2001) 
Caoak et al. (2004^ 
Present work 



-0.45 22.0 < Bj (Vega) < 27.0 

0.442±0.003 19.0 < Bj (Vega) < 23.5 

0.491±0.009 19.0 < Bj (Vega) < 24.4 

0.396±0.001 22.4 < Bj (Vega) < 26.9 

~0.16 26.0 < F450VK(AB) < 29.0 

~0.39 23.0 < F450V1/(AB) < 26.0 

~0.464 16.0 < Bj (Vega) < 21.0 

~0.45 20.0 < fej (Vega) < 25.0 

~0.31 25.0 < F450iy(AB) < 28.0 

0.4±0.1 22.25 < F450W(AB) < 25.25 

0.19±0.01 25.25 < F450Ty(AB) < 27.25 

~0.25 15.0 < B(Vega) < 29.0 

0.473±0.006 16.75 < Bj (Vega) < 24.75 

0.479±0.005 14.4 < Bj (Vega) < 23.6 

0.450±0.008 20.0 < Bj (Vega) < 25.5 

0.497±0.017 21.0 < B(Vega) < 24.5 



Note. — Slopes sivcn in the B-j ba nd have been all taken from 
Table 5 of IKiimmel fc Wagne^ 1200111 . 



6.2. U- and B-band Count Slope 

The slopes of our number count distributions were ob- 
tained by least-squares fits, yielding d\og{N)/dm= 
0.50±0.02 for 5=21-24.5, with x^=0-018, and 
d\og{N)/dm^ 0.48±0.03 for C/=21-24, with x^=0.033. 
In Table [3 our U and B count slopes are listed together 
with those from several surveys. Our slopes are in good 
agreement with othe r studies, with the except ion of 
iWilliams et alJ |)199(tD in both bands and Metcalf e et al.l 
()2001t) in B. The different photometric bands and 
magnitude ranges ca n explain the discrepan cies. More- 
over, differences with IWilliams et al.l l)1996|) could arise 
from the fact that their sample was selected in a red 
band (F606H^-|-F814M^), and hence it is biased towards 
re dder objects; while th e change of the slope at B ^24 
in IMetcalfe et al.l l|200lD fl attens it much more than in 
other studies, as noticed bv lKiimmel fc Wagneil (j2001'). 

Yasuda ct al. (2001) pointed out that, at the bright 
magnitudes where cosmological and evolutionary correc- 
tions are relatively small, the shape of the galaxy num- 
ber counts-magnitude relation is well characterised by 
N{in\) cc 10 the expected relation for a homo- 
geneus galaxy distribution in a Euclidean Universe. But 
since the night sky would be infinitely bright in U and 
B if this trend continued forever, at some f aint magni- 
tude t he U-h&nd count slope must break. iHogg et all 
I^1997^ predicted that [/-band break should happen at 



27 < U < 28, where the median U — R reaches the value 
where o bjects get the UV spec tral slope of a star-forming 
galaxy. iVolonteri et alJ l)2000(l do not find evidence of a 
turn over or fiattening down to F300W^(Vcga)=26, con- 
trary to wh at is claimed by Pozzetti ct al- (,1998). On the 
other hand, Willia ms et alJ l|l99d) report a clear change 
in slope at U ~25.3 on the HDF-N number counts. From 
our data, we report a change of the slope of the U counts 
1.5 mag brighter: U ^ 23.25. Considering that, with our 
combination of depth and area, our survey covers the 
18.25 ^ U < 25.25 mag range with complete data and 
good statistics better than the other existing surveys, 
this brea k should be more sign ificant than the one re- 
ported bv IWilliams et al"! ljl996j) . whose area is less than 
0.5% ours. 

The change of slope at B is reported by ILillv et all 
l)1991l) at B ^ 25. In our covered range, our B num- 
ber counts do not exhibit a clear change. At i? ~ 24, 
it seems that there is a turnover, but we would need 
deeper images in B in order to corroborate it. Nev- 
ertheles s, several studies found this turnover at B ^ 
24 (see iKiimmel fc Wagneil I200H lArnouts et all Il999t 
iWilUams et alJll996HMetcalfe et alJll995(l . so that slope 
change in our B counts is probably real. Perhaps more 
interesting is the fact that from our B data we can con- 
fir m the slight increase in the slope at i? ^ 23, reported 
bv lMetcalfe et all l)1995l) . This couples with the decrease 
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in the slope at fainter magnitudes and leads to an small 
upward bump in the counts centered at B ^ 24, a feature 
that c an be observed in our data. In fact. lMetcalfe et aP 
(|1995j) indicate that this "hump" is a characteristic fea^ 
ture of pure luminosity evolutionary models (PLE) and 
is caused by strongly evolving early-type galaxies at high 
redshift. 

7. MODELING THE GALAXY NUMBER COUNTS 

The large area x depth product of our U and B counts, 
and the simultaneous availability of U, B and K num- 
ber counts for the same field, permit a useful comparison 
with the predictions of galaxy formation models. We re- 
strict this comparison to traditional number count mod- 
els which evolve the z = luminosity function back in 
time, as opposed to SAM models that evolve galaxies for- 
ward in time to z=0. While number count data in blue 
bands can be reproduced with a fairly wide range of num- 
ber count model parameters, reproducing number counts 
in the NIR has proven more challenging, due to the pe- 
culiar change in the number count slope at Kg = 17.5 
(Gardner et al. 1993; CH03). CH03 showed that a fairly 
recent formation epoch for massive galaxies (z ~ 1.5) 
yielded a number count slope change similar to that ob- 
served; allowing red, massive galaxies at higher redshifts 
would yield higher predicted faint counts than seen in 
the data. 

In this paper we extend the number count modeling 
presented in CII03 by comparing the models to both the 
NIR counts and the U and B counts. The combination 
of the blue number count distributions, which are al- 
most featureless over our magnitude range, with the Kg 
number count distribution, which shows a knee at inter- 
mediate magnitudes, provides useful constraints on the 
formation history of the various galaxy types. We pro- 
vide the first number count model that accounts for both 
blue {U, B) and NIR {Kg) number counts. Working on 
the same area of the sky ensures that the distinctions 
between blue and NIR number count profiles are not a 
reflection of cosmic variance. 

To build galaxy numbe r count predicti ons, we have 
used the ncmod code from IGardneil ljl998fl . made avail- 
able by the author at his web site^^. The reader is re- 
ferred to the above publication for details of the model. 
Briefly, the code evolves the local LF back in time, for a 
number of galaxy types, using SEDs from the Galaxy 
Isochrone Synt hesis Spectral Evolution Library (GIS- 
SEL96) model ijBruzual fc CharlotllT993t iLeitherer et all 
Il996j) . The star formation history for each galaxy type is 
parametrized by the redshift of galaxy formation Zf and 
the timescale r of the (exponential) decay of the SFR. 
The code allows for the inclusion of extinction by dust 
internal to the galaxies; dust is modeled as an absorb- 
ing layer, symmetric around the midplane of the galaxy, 
whose thickness is a fraction of the total thickness of the 
stellar disk (Bruzual et al. 1988; Wang 1991). A power- 
law (X extinction law is adopted, and galaxies are 
assumed to have an extinction r45oo = 0.6 [Lz^q/ L*)^-^ 
at 4500 A. (Gardner sets the coefficient in the previ- 
ous equation to 0.2. As discussed below, 0.6 is more 
justified by observations, and we modified the code ac- 
cordingly.) Two recipes are provided to account for the 

http:/ /survey.gsfc.nasa.gov/'^gardner/ncmod/ 



effects of merging, namely, a z-evolution of the LF pa- 
rameters, i.e., $* cx (1 -I- z)^ ^ which conserves the lu- 
minosity density by setting L* oc (1 -I- z )^^, and the 
formulation proposed by Broadh urst et al. (1992), with 
$* (X exp{-(2/^[(l -I- z)"*^ - 1]}. Here, Q is approxi- 
mately the number of objects at z — 1 that will merge 
to form a typical galaxy today, and /3 is a function of the 
look-back time. 

The main inputs of the model are: LFs, SEDs and for- 
mation redshift Zf for each of the galaxy classes, extinc- 
tion and merging switches, and cosmological parameters. 
We have used the local, morpho logically-depende n t lumi - 
nosity functions (MDLF) from INakamura et all lj2003(l . 
which are derived from about 1500 bright galaxies of 
the Sloan Digital Sky Survey (SDSS) northern equatorial 
strips. A literatu re search from 1988 to 2003 showed that 
INakamura et alJ (^003) provide the morphologically- 
dependent LFs with best statistics. After correcting for 
Galactic extinction, the limiting magnitudes are r* < 
15.7 oi B < 17.3 Vega mag. With these depths, the 
MDLFs include local blue compact dwarf galaxies, which 
have typical magnitudes in the rang e B 12—17 
mag. We adopt the galaxy classes from INak amura et alJ 
lj2003l) . who classify galaxies into four groups: E-SO, 
SO/a-Sb, Sbc-Sd, and Im. In Table we give the 
Schechter parameters of these MDLF in the SDSS r* 
filter, for our adopted cosmology. 

We adopt fairly standard population parameters to de- 
scribe each galaxy class (see Tabled. A Salpeter IMF 
is used f or all classes. In view of the mass-metallicity 
relation llTremonti et alJl2004|) . Solar metallicities are 
adopted for the E-SO and SO/a-Sb groups, and lower 
metallicities for later types. Star formation is instan- 
taneous (SSP models) for E-SO, exponentially-decaying 
for spirals, and constant for the Im class. 

We include number evolution, as ample evidence 
shows tha t merger fractions increase with look- 
back time llLe Fevre et all 120001 IConselice et"all I2003t 
iCassata et alJ l2005j) . We parametrize merger-driven 
number evolution as <&* c>c (1-1- z)^ , with /3 — 2.0 
providing good results, and we explore also <&* cx 
exp{— (3//3[(l -I- z)~'^ — 1]}, with Q w 1 giving reason- 
able results^^ 

An accurate description of extinction is critical for 
comparing model predictions with data for wavebands 
from U through Kg. We find Gardner's assumption that 
tb = 0.2 too conservative in view of the evidence from 
studies of local disk galaxies. IPeletier fc W illne? (!1993) 
infer a face-on tb ~ 1-0 for spirals: [ Boselli fc Gavaz^ 
(199?) find TB > 1, while iXilouris et all l| 19991 infer a 
central tb = 0.8. We have adopted tb = 0.6, an interme- 
diate value between Gardner's and the above references. 
Extinction is assumed for all galaxy classes. 

7.1. Number Count Model Predictions 

When evolving galaxy populations back in time, a 
characteristic power-law behavior is found for the counts. 
This is given by the intrinsic brightening of galaxies 
with look-back time together with the z-evolution of 
the volume element, and it extends faintward until we 

Note that the choice of exponents is not based on exponents r) 
found when parametrizing merger fractions as fm ~ (1 -I- z)^'- the 
latter does not lead to a ^-dependency of <I>* of the form oc (1 + ^)'^. 
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Fig. 10. — From top to bottom, evolution models overplotted 
on the {/, B, and galaxy number counts, with Zf = 1.5 for 
ellipticals and early spirals and Zf = 4 for other types. Number 
evolution of the luminosity function is modeled using <i>* oc (1+z)'', 
with /3 = 2.0. Fill ed c ircle s: Our GWS data. U and B counts 
are from Tables 1121 and IlijI respectively, and Kb counts are from 
CH03. Open circles: Counts from other authors. See references in 
the text. Solid lines: Total galaxy count prediction. Dotted lines: 
Prediction for ellipticals. Short-dashed lines: Prediction for early 
spirals. Long-dashed lines: prediction for late spirals. Long-short 
dashed lines: Prediction for Im galaxies. 
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Fig. 11. — From top to bottom, evolution models overplotted 
on the U, B, and Kb galaxy number counts, with Z{ = 1.5 for 
ellipticals, and 2f = 4 for other types inc luding early- type s pirals . 
Number evolution is modeled using the I Broadhurst et all <1992l) 
prescription, with Qj, = 1.0. Fille d ci rcles : Our GWS data. U 
and B counts are from Tables 1121 and IlijI respectively, and Kb 
counts are from CH03. Open circles: Counts from other authors. 
See references in the text. Solid lines: Total galaxy count predic- 
tion. Dotted lines: Prediction for ellipticals. Short-dashed lines: 
Prediction for early spirals. Long-dashed lines: prediction for late 
spirals. Long-short dashed lines: Prediction for Im galaxies. 
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TABLE 8 

SCHECHTER PARAMETERS OF THE ADOPTED LUMINOSITY FUNCTIONS 



LF Parameters M*(r*)''''' a #* (0.001 Mpc'^) b 



(1) (2) (3) (4) 



E-SO -21.53± 0.17 -0.83±0.26 1.61 ± 0.09 

SO/a-Sb -21.08 ± 0.19 -1.15±0.26 3.26 ± 0.15 

Sbc-Sd -21.08± 0.20 -0.71±0.26 1.48 ± 0.05 

Im<^ -20.78 -1.90 0.37 



''M*(7-*) in AB magnitudes IStoughton et aljl2n()3l . 
''Converted to ?i = ffo/lOO = 0.7. 
'^Irregular, star-forming galaxies. 



TABLE 9 

Adopted Properties for each Galaxy Type 





SFR 








Galaxy Type 


Type 




Z/Z0 


IMF 


(1) 


(2) 


(3) 


(4) 


(5) 


E-SO 


SSP 





1 


Salpeter 


SO/a-Sb 


exponential 


4 


1 


Salpeter 


Sbc-Sd 


exponential 


7 


2/5 


Salpeter 


Im 


constant 




1/5 


Salpeter 



^ The star-formation c-folding time in Gyr. 



reach magnitudes for which galaxies have yet to build a 
large fraction of their final stellar masses. Breaking such 
power-law behavior, as is necessary at intermediate mag- 
nitudes for the Ks counts, can only be accomplished by 
setting the formation redshift zt of a dominant popula- 
tion to a moderately low value; other model parameters, 
such as the evolution of the merger fraction or the ex- 
tinction, are unable to yield the observed knee in the 
Ks counts. A model that reproduces the [/, B, and 
Ks counts can be obtained by setting Zf = 1.5 for el- 
lipticals and early spirals, and Zf = 4 for other galaxy 
classes, and assuming a merger-driven number evolution 
as (1 -|- z)^'*^'. Figure [Till shows the number count pre- 
dictions of this model vaU , B and Kg. Our total num- 
ber counts as well as the contributions of each galaxy 
class are shown. Number count data from other au- 
thors have been plotted: in the J7-band fVolonter i et al.l 
EpOO; Metcalfe et al. 2001; Yasuda et al. 200 j); in 
the B-band ( Gardner et al.j|199^yyolonteriet_ayj200^ 

[Metcalfe et al 2001.: .Yasuda et all I200H iHuang et all 

I2001bt IKiimmel fc Wagned 120011): and in the K,- 
band ("Gardner et al.' '1993"; 'Diorgovski et al." 119951 
[McLeod et al. 199 5; Hu ang et al. 1997; Minczaki et all 
IT998HTotani et al.ll200l() . 

A second model which avoids the uncomfortably low Zf 
for spirals is shown in Figure^] This model, with Zf = 4 
for early-type spirals, also provides a good description 
of the observe d counts in U , B^ and i vTs. The merger 
prescription bv iBroadhurst et alJ 1)199 2') was used, with 
Q = 1, in order to lower the slope in the faint Ks counts, 
which seems steeper than the trend defined by the data; 
shallower faint-end slopes are obtained for lower values 
of Q (§01. For ellipticals, we need to keep Zf = 1.5, oth- 
erwise we fail to reproduce the — 17.5 knee (S I7.1.2() . 
While differing in the details, both models presented in 
Figures E| and ^] provide similarly accurate overall de- 



scriptions of the data. 

These models capture the essential ingredients of the 
modeling of counts we presented in CH03, by setting 
Zf to a low value, and indicate that the solution proposed 
by CH03 to the knee in the A's counts docs reproduce the 
counts in blue passbands as well. But the new models 
differ from CH03 in several details, including the LF, 
here derived from SDSS data; the addition of extinction, 
which was unnecessary for the NIR counts of CH03; the 
inclusion of merging evolution; and the exclusion of any 
galaxy population not present in the local LF. 

7.1.1. Dominant populations 

Over the magnitude range sampled by our data, the 
UB counts are dominated by early- and late- type spi- 
rals, especially so in U . At bright magnitudes, the lower 
contributions of E-SO and Im galaxies reflect their lower 
4>* at z = (see Table|Hl). Im galaxies contribute ^^1/20 
of the total counts throughout our magnitude range, 
while the contributions of E-SOs decreases with mag- 
nitude owing to the strong K-corrections of their red 
stellar populations: neither has a measurable contribu- 
tion to the counts. In A's, the counts are dominated 
by E-SO at bright magnitudes, and by early-type spirals 
at Ks > 19.5. We note that, c ontrary to many previ- 
ous works ( fBabul fc Reeij'1992': 'Philhos fc Driver''l99| 
Driver et al . 1995; Babul fc Ferguson 1996; EUis 199^ 
Driver et al . 1998; Phill ips et aL 2000), we do not need 
to introduce ad-hoc populations of star-forming galaxies 
at intermedite redshifts to explain the observed number 
counts; our model counts only comprise galaxies in the 
local LFs evolved back in time. 

7.1.2. Formation redshifts 

As discussed in CII03, only by adopting a late Zf ~ 1.5 
for clhpticals do the models reproduce the Ks = 17.5 
knee. In the context of the number count models we 
are using, this feature requires the disappearance, or sig- 
nificant downsizing, of a dominant red population with 
look-back time at intermediate redshifts. Only the el- 
lipticals can play this role given the z = LFs. We 
searched for alternative schemes which did not involve a 
late Zf for ellipticals, by varying either the star forma- 
tion time scales, the merger recipes, or the extinction. 
We found no combination of parameters that simultane- 
ously yields the knee in Ks and the lack of a knee in blue 
counts. In particular, postulating a specific dusty phase 
in galaxy evolution at intermediate redshifts barely af- 
fects the Ks counts, and instead introduces features in 
the blue number counts which are not seen in the data. 

The merger rate in hierarchical ACDM models, and the 
merger-driven star formation assumed to derive f rom it 
in SAM models, peaks at z ~ 2 (e.g.. lCole et al.i r2000). 
Ellipticals are assumed to result from such mergers. We 
conjecture that the knee at Kg = 17.5 in the NIR num- 
ber counts may reflect the onset of this red population 
of merger origin. Such merger origin cannot be properly 
accounted for by ncmod, due to its simplified galaxy for- 
mation recipes, in which mergers change galaxy numbers 
and luminosities, but not galaxy types. The appearance 
of the red population needs to be introduced by setting 
the appropriate Zf. Setting the population to form in- 
stantaneously (SSP model) helps ncmod mimic the onset 
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of the merger-induced red population by minimizing the 
phase during which the galaxies contain young stars. 

The late (z ^ 1.5) formation of ellipticals predicted 
by ncmod can be reconciled with the existence of a sig- 
nificant po pulation of evolved galax ies at 2 < z < 4 
llFranx et al. 2(m: T.a,hhI7ra,].li2nn,'^ IDaddi et HnM)^ 
iGlazebrook et al.ll2004Hde Mello et al.ll2004|) . The mass 
density of red, dead galaxies at z=1.8 was ~16% of the lo- 
cal density (McCarthy ct al. 2004); Zf — 1.5 in our num- 
ber count models reflects such downsizing of the evolved 
population. Also, the rapid increase in the clu stering of 
high-z galaxies toward redder colors ijDaddi e t al. 2003) 
suggests that red sources at 2 < z < 4 trace a clus- 
ter population; our Groth number counts predominantly 
trace a field population, in which elliptical formation 
may occur later than in clusters. A similarly late forma- 
tion for the majority of ellipticals/SO (zf < 2) is derived 
from the analysis of color gradients in early-type galax- 
ies in the HDFN using m ulti-zone single collapse models 
ijMenanteau et al.ll200ll) . 

In contrast with the ellipticals, Zf for all the star- 
forming galaxy classes is poorly constrained in the mod- 
els. For early spirals, values in the range 1.5 < Zf < 10 
yield accepable results; for late spirals, the acceptable 
range is 0.5 < Zf < 10. 

7.1.3. Dust 

Adoption of extinction for all galaxy classes is criti- 
cal to the success of the model. In particular, we need 
to assume extinction for ellipticals. Otherwise, the pre- 
dicted blue counts show a strong feature, at [/ ~ 18 and 
B ~ 19. Our extinctions {tb — 0.6 for L* galaxies, see 
§[7|l are moderate, if we take into account that precur- 
sors of ellipticals may include dusty EROs at 1 < z < 2 
(e.g. lSmail et aljr200 ^. Red, massive galaxies at higher 
redshifts show much higher exti nctions; e.g.. Ay ^ 2.7 
at 2 < z < 3.5 in the HDF-S ijForster Schreiber et all 
l2?ffll . 

7.1.4. Color distributions 

As pointed out by iGardneJ l)1998D , color distributions 
provide checks on the evolutionary processes of galaxies 
that cannot be obtained from number counts. We fo- 
cus on the V ~ I color, as our observed catalog is most 
complete in those bands. The analysis starts with our 
model with Zf = 1.5 for spirals, and focusses on galaxies 
with 21 < y < 23; similar conclusions are derived from 
other magnitude ranges. V — I color distributions for the 
model and the data are compared in FigurelT^. Clearly, 
the model distribution is narrower than seen in the data, 
and it shows a deficit of red galaxies. Photometric er- 
rors cannot account for the discre pancy, as s hown by the 
dotted distribution in Figure di. iGardneil ((1991 noted 
that the color predictions from the models are often too 
narrow owing to the discrete nature of the modeled pop- 
ulations. We find that this is the case for spirals (see 
Fig. 112b. where we have plotted early- and late-type spi- 
rals together) and for irregulars fFig. 112111. However, the 
model distribution for E-SO is broad (Fig. Eb)- Such 
width arises from the strong K-corrcctions of old popu- 
lations, which appears to be responsible for the width in 
the observed V — I histogram. 

Figure El suggests that the discrepancy between ob- 
served and model colors arises from the failure to include 



red, evolved populations in the precursors of present-day 
early-to-intermediate type spirals. Possibly, our early- 
spiral class (SO/a - Sb) is too broad to account for the 
range of colors present in galaxies of types SO/a to Sb. 
The models appear to miss the strong contribution of 
bulges to the integrated light of early-type spirals. We 
are constrained by the morphological classes of the lu- 
minosity function defined by Nakamura et al., hence we 
cannot properly split the early S class. To estimate the 
potential effects of splitting the early S class, we generate 
a toy-model in which we simply divide the early spirals 
in two equal groups. For the first ("SO/a-Sa" group) we 
assign the formation history corresponding to ellipticals 
(SSP Solar mctallicity, Scalo IMF, and formation redshift 
Zf = 1.5). For the second ("Sb group"), we assign a for- 
mation redshift Zf = 4.0, an exponentially-decaying SFH 
with T = 7 Gyr, and Solar metallicity. The color distribu- 
tion for such model (Figure E^), provides a close match 
to the observed colors, specially when smoothed to ac- 
count for photometric errors (dotted line in Figure Et)- 
This exercise suggests that the deficit of red galaxies in 
the model shown in Figure 112b traces an overly simpli- 
fied description of the populations in early-type spirals. 
The single-age nature of the population might affect as 
well; a broad color distribution for the S group would 
be obtained as well if an important fraction of the pre- 
cursors of present-day spirals became dominated by old 
populations earlier than assumed in our models. Finally, 
we note that splitting the early-S population does not 
affect significantly to the fits to the number counts. 

7.1.5. Cosmological parameters 

When the cosmology is changed to an Einstein-de Sit- 
ter model, our ncmod model reproduces the U and B 
counts, but fails to reproduce the Kg counts by a large 
margin. Hence, the combined number counts in UBK 
appear to favor a A-cosmol ogy over an Einstein-de Sitter 
one, a result also found bv lTotani fc Yoshiil l|200(H) . 

8. SUMMARY 

We have presented U and B number counts from a field 
of the GOYA Survey that covers ^--^900 arcmin^ over the 
Groth- Westphal Strip. Counts are derived from 0.2350 
deg^ and 0.2467 deg^ of sky in U and B, respectively. 
Achieved limiting magnitudes (50% detection efficiency 
for point sources) are U — 24.8 mag and B — 25.5 mag, 
in the Vega system. The counts have been corrected for 
detection efficiency as a function of source size, and for 
spurious detections, using the method of S/N thresh- 
old in two complementary half-time images detailed in 
Elichc-Moral & BalccUs (2005. hereafter EB05). Star- 
galaxy separation has been performed using HST images 
in the GWS and stellarity indexes from SExtractor. 

Counts are given over 18.0 < U < 25.0 and 19.5 < 
B < 25.5. These wide ranges (7 mag in U and 6 mag in 
B) result from the combination of wide area and depth 
of our survey. In both bands, our number counts are in 
good agreement with other studies that cover fainter and 
brighter magnitudes. The slopes of the number count 
distributions are very similar to those reported by other 
authors: d\og{N)/dm = 0.50 ± 0.02 for 5=21.0-24.5, 
and d\og{N)/dm = 0.48 ± 0.03 for [7=21.0-24.0. 

When combined with number counts, the data pro- 
vide strong constraints on galaxy formation models, due 
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Fig. 12. — The V — I color distributions for the data {hashed histogram), overlaid on four representations of the color distributions for 
the model with Z{ = 1.5 for spirals {Thick-lined empty histogram). The model is normalized to the number of galaxies in e3x:h V magnitude 
bin. Only galaxies with 21 < y < 23 are shown. Panel (a): color distribution for the entire 21 < V < 23 model population. The effect of 
photometric errors has been modeled by convolving the model histogram with a Gaussian kernel of ayi = 0.13 {dotted line). Panel (b): 
Color distribution for E/SO. Panel (c): Color distribution for SO/a-Sd. Panel (d): Color distribution for Im. 



to the presence of a knee at Kg = 17.5 in the NIR 
counts, and the absence of such feature in blue pass- 
bands. Adopting a A-dominatcd cosmological model 
(flM = 0.3, flA ^ 0.7, Ho = 70 km s"! Mpc"!), a simple 
number count model including luminosity evolution from 
GISSEL SEDs and number evolution as (1 + 2;)^-^ accu- 
rately reproduces the observed counts in f/, B, and 
in a consistent way. Extensive modeling suggests that, 
only by assuming a moderately low formation redshift 
(zf « 1.5) for the dominant NIR population (ellipticals) 
does the model reproduce the = 17.5 knee; while, 
reproducing the lack of a knee in U and B counts in 
turn requires the adoption of a moderate optical depth 
for all galaxy types, including ellipticals {tb = 0.6 for L* 
galaxies). Neither of the two assumptions is at odds with 
current ideas on galax;y formation and evolution in hier- 
archical Universes. Future comparison with SAM galaxy 
formation models will tell whether the Zf — 1.5 required 
by the model reflects, as we suspect, a major epoch of 
early-type field galaxy formation through mergers of disk 
galaxies, rather than the epoch for the formation of the 
stellar content of ellipticals. 
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Fig. 13. — The V — I color distributions for the data {hashed histogram), overlaid on four representations of the color distributions for the 
toy model in which the early-type spiral class has been split in two halves, corresponding to ad-hoc SO/a-Sa and Sb populations, respectively. 
The model is normalized to the number of galaxies in each V magnitude bin. Only galaxies with 21 < V < 23 are shown. Thin-lined filled 
histogram: Color distributions of our data. Thick-lined empty histogram: Color distributions predicted by the model. Panel (a): color 
distribution for the entire 21 < V < 23 model population. The effect of photometric errors has been modeled by convolving the model 
histogram with a Gaussian kernel of ayi = 0.13 {dotted line). Panel (b): Color distributions for the ad-hoc SO/a-Sa population. Panel 
(c): Color distributions for the ad-hoc Sb population. Panel (d): Color distributions for Sbc-Sd. For this model, the color distribution for 
E-SO galaxies is identical to that shown in Fig. 112b . 



APPENDIX 

EFFICIENCY ERROR ESTIMATION BY FUNCTIONAL AND MATRICIAL METHODS 

Our error treatment is analogous to that described in the Appendix of CH03. Number count errors arise from 
detection efficiency errors and statisti cal counting, once the spurious have been removed from the catalogues. The 
latter errors can be approximated by ijGehreyil 986(1 : 



a^f ^*(mout) = ■\/iVdct,.(™out) + 7 (Al) 



4 



crf/*''*(rnout) = Adot,s(mout) - A"dct,s (^out) 



H 3 



9^det,s(TOout) 3v/Afdct,s(mout) 



with s being the index that refers to the size group (see ti5.1|) . Adct,s('7iout) the raw number counts in mout for the 
size group s after spurious subtraction, ct^ the upper error, and crj^j''*^* the lower error. Depending on the method 
used for completeness correction, errors have been estimated as follows: 

1. Using the method of functional efficiency, E{mout)- Considering the counting errors from equations (|A1|) and 
erro rs arising from t he efficiency corrections, la error propagation of equation Q for the size group s leads to 
(see lBevingtonlll969D 
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(A2) 



where ^^*^'(mout) is the statistical error associated to counting (upper or lower), £'(mout) is the efficiency at 
that magnitude bin for the size group s, A is the area of counting, and cr^(TOout) is the efficiency error. The 
factor 0.5 accounts for the magnitude interval used for counting (0.5 mag in our case). 

Using upper and lower statistical errors from equations HA1|) . we can obtain the corresponding total upper and 
lower errors af^{mout) and erf j (wout ) ■ Then the propagated errors from the addition of the number counts of 
the three size groups are 



O'l.totall'Tlout) — 



Cru,total(^out) — 



\ s=l 



(A3) 



Finally, this error is quadratically added to that associated to star counting (see the definition of this last in 



CrLfinaKmout) = V [cri,total(TOout)]^ + [cri,stars(TOout)]^ 



(A4) 



Cru,final("T.out) = \/[cru,total("T-out)]^ + [Cu.stars (Wout )] " 



2. For the method of the efficiency matrices, P(min;TOout)- We have considered counting and efficiency errors 
separately: 

• We have estimated statistical counting errors by adding lower and upper statistical errors from equation 
iPTjl to the raw counts in equation (|10|l in each size group: 



^^u,s(min)' 




[iVdot,s(TOout) + crf;f''*(mout)] Ps ^(mout,min) 



(A5) 



A^dot,s(TOout) + Crf/*'''(TOout) • -Ps ^("^out, TOin) 



0.5 • 



where iVi,s(win)' and Afu,s(win)' sue the lower and upper values to the corrected number counts in equa- 
tion H10|l respectively; -/Vdct,s(™out) is the raw number counts in mout after spurious subtraction; and 
P~^(mout, 'Tiin) is the element (mout, Win) of the inverse of the efficiency matrix that corresponds to the 
size group s, Ps(min; mout)- Notice that the addition of the upper limit to the raw counts would produce 
the lower limit in the efficiency-corrected counts, because of the inversion of the matrix, and vice-versa. 
Therefore, upper and lower errors are defined as the difference of the resulting values with respect to those 
obtained without adding errors: 



Crfs(TOin)' = Norig.simin)' - NisimnY 
O'us('^in)' = A^u,s(min)' - Aforig,s(TOin)', 



(A6) 



where A'orig.s(min)' is obtained from equation H10() . 

The errors from efficiency correction are estimated in a similar way: upper and lower limits to the efficiency- 
corrected counts arc computed by adding or subtracting the RMS of the efficiency matrices for the size group 
s in equation pO|l . Then, corresponding upper and lower errors due to efficiency determination errors are 
the difference of these resulting values with respect to those obtained in equation Hl()|l . as follows: 



^^I("^in)' = ^orig,s(TOin)' " qXA ' {X]vmout -^dct,s (^out ) ' (Ps + ^s) H™out,mi„)j | 
a^^^^min)' = ■ {5I)vm„„t A^dct,s(mout) • (Ps - ^p7s)"H'^out,min) | - A'^orig^s (win)', 

where ixp^ is the RMS matrix of the efficiency matrix for the size group s (see Figure EJ. 



(A7) 
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TABLE BIO 

Relative Positions and Rotations of INT/WFC CCDs 



CCD ID 


Ori 


jinal coefT. 




Transformed coefT. 


Corrected coefF. 


Offsets 














x^ 


Y^ 


AXi iXYi 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) (10) 


#1 


0.01868 


-336.74 


3039.14 


4338.73 


61.97 


4339.73 


57.47 


1.0 -4.5 


#2 


-90.62115 


3180.68 


1729.67 


2305.90 


6298.23 


2308.90 


6296.73 


3.0 -1.5 


#3 


0.11436 


3876.73 


2996.30 


130.26 


96.97 


134.76 


99.47 


4.5 2.5 


#4 


0.00000 


1778.00 


3029.00 


2223.00 


72.00 


2223.00 


72.00 


0.0 0.0 



Note. — Col. (2)-(4) a re or iginal values for the rotation angles in degrees and positions of the optical axis in the frame 
of each CCD, reported bv lTav lor ( 2Q0Q^. Col. (5)-(6) arc these positions in a global reference frame where the optical axis 
is situated at coordinates (2223.0,72.0). Col. {7)-(S) are our corrected coefficients in the same global frame. Col. (9)-{10) 
show the computed offsets. All positions are given in pixels. 

Thus, the errors from equations ljA6p and (|A7p are added quadraticahy in each case to obtain the upper and 
lower values: 



Cu^sC^JT-in) = \J [o-J^,s(™in)']^ + (»^in)' 

o-i,s(™in) = \ o-pf(min)' + (rf^imi^y 



2 

I (A8) 

The errors associated to star counts are estimated in an analogous way. The quadratic addition of values from 
this last expression for the three size gro ups a nd th e errors associated to star count subtraction would estimate 
the final total error, as in equations HA3() and HA4|I . 

INT/WFC OFFSETS 

In ti4.3l we report the difficulties we found when we tried to fit an independent astrometric solution for each CCD 
of the INT/WFC, using the TNX projection from IRAF astrometric tasks. The instrument exhibits an important 
"pincushion" distortion introduced by the telescope optics (see an instrument description in SQ- We concluded that 
a single fit to the entire field of view was needed in order to produce a continuous solution in t he CCD edges. 

Our procedure for creating single images of each exposure is based on the method described bv lTavloJ ( 2000). He did 
an astrometric solution of the INT/WFC, correct to an accuracy of 1 or 2 pixels 0.5"). The corrected coordinates 
(a;/, yt) were obtained from each set of pixel coordinates (x^, ?/j) by first translating so that the origin is on the optical 
axis, then rotating to the correct angle, and finally correcting for the radial distortion effect, as follows: 



(Bl) 



where {Xi, Yi) are the coordinates of the optical centre of the instrument in the pixel coordinate system of CCD#i, 
9i is the angle at which CCD#i sits on the focal plane, and D is the pincushion distortion c oeffic ient. The values of 
the coefficients in these equations reported by Taylor (2000) are listed in Col. (2)-(4) of Table IbToI The optical axis is 
the final origin of the reference frame. 

Single images of each exposure were created by posi tioning each C CD in an empty frame wide enough to include 
the four chips, and using the coefficients reported bv iTavlon l)2000() . The optical axis was situated in coordinates 
(2223.0,72.0) in the global frame, and coordinates (X^, Yi) were transformed to this new origin, ta king i nto account 
the angles between the CCDs. Taylor's transformed coefficients are tabulated in Col. (5)-(6) of Table IbTIiI After that, 
we fitted a third-order astrometric solution with the TNX projection, in order to describe remaining, non-corrected 
linear distortions and non-linear optical aberrations. Disconti nuities in the astrometric solutions indicated that errors 
of ^l"are present in the chip separations reported bv iTavloJ (|2000). 

In order to improve the astrometric solution, we estimated corrections to Taylor's positions by offseting the relative 
CCD positions until we found a combination which minimized the astrometric discontinuities at the chip-chip interface. 
The corrected values and offsets are listed at Table IbToI 

Discontinuities in the astrometric solutions were removed, and the RMS astrometric error after a new third-order 
fit to the astrometry has been reduced to 0.3"over the ^36'x36' field of the camera. 
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TABLE 11 

Estimated Fractions of Real Sources Rejected as "False Detections", Truly Spurious Sources Rejected as "False Detections", and Truly Spurious 

Non-Rejected from the U Catalogue, for Different DETECT_THRESH-(5/W)iini Combinations 
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Note. — All percentages are relative to the total number of detections in the U catalogue. The fraction of non-rejected truly spurious sources (last row) has been computed subtracting the number of rejected truly- 
spurious detections (penultimate row) from the estimated maximum number of truly spurious sources that are in the catalogue, which is --^ 29% of the total number of detections for DETECT_THRESH= . 4 , 0.5, and 0.6, and 
25% for DETECT_THRESH^0.7 (see the text). 



TABLE 12 

U Differential Number Counts in the GWS Field 
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Note. — Col. (2), (3), and (4) show source numbers (raw counts) for size groups rc < 1.5 , 1.5 < < S", and ^e > 3 , respectively, in intervals of 0.5 mag centered on the 
m.agnitudes shown in Col. (1). Col. (5), (6), and (7) show raw counts corrected for spurious detections. Col. (8), (9), and (10) indicate the applied functional efficiency correction for 
the spurious-corrected number counts in each one of the previous size groups. Col. (11) presents differential galaxy number counts per unit magnitude and area (spurious-corrected, 
efficiency-corrected, and with stars subtracted). Col. (12) and (13) show the 1-a confidence upper and lower limits on corrected differential counts of Col. (11). Col. (14) lists the 
logarithm of Col. (11). 



TABLE 13 

B Differential Number Counts in the GWS Field 



to 
O 



B 

(1) 


Rawl 
(N) 

(2) 


Raw2 
(N) 
(3) 


Raw3 
(N) 
(4) 


(nT 

yiy J 
C5) 


Ar2 
(6) 


Ar3 

spur 

(TV) 
(7) 


Ei{U) 
(8) 


E2{U) 
C9) 


fio^ 


N 


0"u 

f TV ma[i~^ deff""^) 
("121 


(7\ 

( M iTiae'""'^ dpe"""^^ 
(13) 


log(N) 
(M) 


19.75 


64 








64 










1.000 


1.000 


1.000 


206.90 


238.93 


178.96 


2.316 


20.25 


82 








82 










0.988 


0.979 


0.970 


480.55 


186.23 


124.56 


2.682 


20.75 


130 


1 





130 


1 







0.981 


0.969 


0.962 


689.34 


233.03 


173.83 


2.838 


zl.zo 


zuo 




U 











U.y / o 


u.yoi 


U.y41 


yoo.oo 




zo / .uy 




21.75 


302 


5 





302 


5 







0.961 


0.936 


0.927 


1902.86 


308.81 


253.34 


3.279 


22.25 


484 


16 





484 


16 







0.948 


0.909 


0.899 


3563.79 


353.15 


301.57 


3.552 


22.75 


808 


70 





808 


70 







0.934 


0.891 


0.878 


6883.04 


439.10 


393.41 


3.838 


23.25 


1267 


213 





1267 


212 







0.922 


0.865 


0.847 


12073.11 


583.52 


542.64 


4.082 


23.75 


1941 


581 


3 


1941 


581 




3 


0.906 


0.841 


0.808 


22130.66 


864.11 


833.89 


4.345 


24.25 


2741 


1114 


4 


2741 


1114 




3 


0.884 


0.812 


0.744 


34401.66 


1258.68 


1230.61 


4.537 


24.75 


3377 


1747 


5 


3376 


1743 




4 


0.856 


0.677 


0.483 


50142.56 


2040.21 


2017.83 


4.700 


25.25 


3891 


2175 


4 


3884 


2140 




4 


0.779 


0.319 


0.211 


92527.11 


7562.70 


7551.66 


4.966 


25.75 


3520 


2382 


11 


3505 


2134 




9 


0.474 


0.164 


0.136 


163445.76 


26245.47 


26235.68 


5.213 



Note. — Columns arc as in Tabic [T2l 



